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1. INTRODUCTION 

As part of the Space Shuttle Return to Flight (RTF) program, Marshall Space 
Flight Center (MSFC) performed computational fluid dynamics (CFD) analysis to define 
the velocity flowfields around the Shuttle stack at liftoff. These CFD predicted velocity 
flowfields were used in debris transport analysis (DTA). High speed flows such as 
plumes induce or ‘entrain’ mass from the surrounding environment. Previous work had 
shown that CFD analysis over-predicts plume induced flows. Therefore, the DTA would 
tend to 1) predict more debris impacts, and 2) the debris velocity (and kinetic energy) of 
those impacts would be too high. 

At a November, 2004 peer-review it was recommended that the Liftoff DTA team 
quantify the uncertainty in the DTA caused by the CFD’s over prediction of plume 
induced flow. To do so, the Liftoff DTA team needed benchmark quality data for plume 
induced flow to quantify the CFD accuracy and its effect on the DTA. 

MSFC’s Nozzle Test Facility (NTF) conducted the “Nozzle Induced Flows test, 
P#2456” to obtain experimental data for plume induced flows for nozzle flow exhausting 
into q quiescent freestream. Planning for the test began in December, 2004 and the 
experimental data was obtained in February and March of 2005. The funding for this test 
was provided by MSFC’s Space Shuttle Propulsion Systems Integration and Engineering 
office. 


2. TEST PROGRAM 

The NTF normally tests subscale nozzles in a sealed test cell that is evacuated to 
simulate altitude. The test medium, or nozzle’s working fluid, is warm air. The NTF was 
modified to move the nozzle test articles downstream and outside of the test cabin to 
enable measuring the entrainment of the surrounding air. See Section 2.2 for a 
description of the NTF and modifications to it for this test. 

This experiment was run with two existing test articles; a Stratford nozzle and a 
Bell nozzle. The Stratford nozzle has a design area ratio of 1.115; the Bell, 12.2. These 
test articles are fully described in Section 2.1. 

The Bell nozzle entrainment data was of primary interest to the Liftoff DTA team. 
Within the limits of the NTF capabilities, it provided a plume structure most similar to a 
Space Shuttle main engine’s plume at sea level. The Stratford nozzle, being close to a 
sonic nozzle, provided an entrainment dataset with simpler plume physics. 

The test data acquired and documented in this report consists of; Schlieren images 
of the plume structure, fog based flow visualization, nozzle wall pressures (Bell only), 
particle image velocimetry (PIV), and hot film anemometry (HFA). The test data 
systems and test setup are described in Section 2.3. The test results are presented and 
discussed in Section 3. 
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2.1. Description of Nozzle Test Articles 

The nozzle test articles are shown in Figure 1. Both are aluminum. The Stratford 
nozzle cross-section is shown in Figure 2. It has a design area ratio of 1.115. Its as- 
designed throat diameter is 35.56mm. The Stratford nozzle’s throat was measured at 
35.725mm. The difference between the as-designed and the as-built throat diameters 
results in a 1% increase in throat area. The nozzle’s internal contour shown on the 
drawing, the as-designed contour, (Figure 2) was generated with the formula as defined 
in Reference 1. Appendix 1 contains a table (Table Al) of nozzle wall points calculated 
with the formula in Reference 1. The Excel® file (‘NozzleIntemalGeo.xls’) used to 
generate the nozzle internal contour in Table Al and for the test related computational 
fluid dynamics is included in the electronic media associated with this report. 

The Stratford nozzle did not have internal wall pressure measurements. 

The Bell nozzle cross-section is shown in Figure 3. It has an area ratio of 12.2. 

Its throat diameter is 39.218mm. The nozzle wall is 12.7mm thick. The as-designed 
internal contour points for this nozzle are in Table A2 in Appendix 1. This internal 
contour is contained within the Excel file, ‘NozzleIntemalGeo.xls’ , which is included in 
the electronic media associated with the report. 

The Bell nozzle had about 40 static pressure measurements on its internal wall. 
The locations of these are also defined in Table A3 in Appendix A. 



Figure 1. Nozzle Test Articles. 
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Figure 3. Bell Nozzle Cross-Section. 
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2.2. Test Facility 

The NTF is designed to determine nozzle performance across a wide range of 
pressure ratios. The supply system can provide up to 5.4 kgm/s of heated (up to 422°K) 
air at 20.4 atmospheres. The test cabin pressure can be reduced to pressures as low as 
6.8e-3atm and is controlled with a two-stage air ejector system. Thrust and off-axis loads 
generated by nozzles are measured with flexure-isolated load cells. Induced loads from 
thermal growth and bringing flow across the metric plane are accounted for with a tare 
process. 

Facility instrumentation used in this test consisted of mass flow and nozzle inlet 
total temperature and pressure. Schlieren images of the plume structure were recorded on 
a video camera. 

The data acquisition system is composed of thermocouples and individual 
pressure transducers from various manufacturers, two Pressure Systems, Inc®. (PSI) 
Model 780B pressure scanners, a Hewlett-Packard® Model 3852 Integrating Voltmeter, 
and two National Instruments® SCXI multiplexers connected to National Instruments 
PCI-MIO-16 cards. Data is collected by three desktop computers. 

2.2.1. Modifications for Current Test 

The NTF was modified, as shown schematically in Figure 4, to move the nozzle 
out of the test cabin to a new test position downstream and outside of the test cabin. 
Photographs of the NTF before and after the modification are shown in Figure 5 and 
Figure 6. To move the nozzle downstream, a pipe was fabricated and attached to the flow 
straightening section that remained in the test cabin. The load cell normally used to 
measure thrust remained in place, but thrust was not a desired measurement in this 
experiment. Furthermore, securing the extension pipe in place invalidated the thrust 
measurement. No valid thrust measurements were obtained. 

The ejector pipe was cut as shown in Figure 6 to become the ‘plume capture 
pipe’. The ejectors were not run during this experiment. However, the ejector hardware 
remained within the plume capture pipe. This hardware is approximately three meters 
downstream of the pipe entrance. 
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Figure 4. Schematic of Test Set Up. 
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Figure 5. Nozzle Test Article in Normal Test Position. 
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Figure 6. Stratford Nozzle in New Test Position. 
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2.2.2. Geometry of Current Test Configuration 

The layout of the room (A-bay of building 4777) containing the NTF is shown in 
Figure 7. There were several small holes in the walls, and three large openings in the 
building though which outside air entered the building during the experiment. The large 
openings were the roof vents, the rollup garage door and the doorway between A-bay and 
B-bay. 


There was a fourth large opening for part of the testing. A window 
(approximately 813mm wide by 508mm tall) in the south wall of the A-bay was open 
during acquisition of the Stratford PIV windows A, B, E, F and C. Between the window 
and the test cabin, a screening panel was placed to divert the inflow and reduce its 
northerly momentum. After the Stratford PIV C window data was obtained, the window 
was covered with a tarp, which significantly reduced the air flow through it. Partially 
covering this window with the tarp was one of the parametrics done to improve the seed 
quality for the PIV. The tarp remained in place through the rest of the PIV data 
acquisition and all of the HFA data acquisition. 

The ceiling vents were approximately one third of the length of the building, 
centered, north to south, in the ceiling. They were about 300mm wide, east to west and 
were opened, vertically, approximately 200mm. The rollup door was 4.2m wide and was 
raised approximately 700mm. The doorway between the A and B-bay is 2108mm tall by 
84mm wide. 

Figure 7 indicates the plume capture pipe centerline was on the same centerline as 
the test articles. However, the plume capture pipe centerline was offset approximately 
25mm south of the nozzles’ centerlines. The plume capture pipe centerline was 
approximately parallel to the nozzles’ centerlines. The centerlines were coincident in the 
horizontal plane (i.e., relative to the floor). 

The Bell nozzle setup is shown in Figure 8 and Figure 9. The Stratford nozzle 
setup is shown in Figure 10. Figure 8 and Figure 10 are both looking south. Those 
dimensions not shown in the Stratford image can be found in the Bell nozzle images. 

The distances from the floor to the outer diameter of the nozzle (in the 1 80° plane) were 
1121 and 1 182mm for the Bell and Stratford nozzles, respectively. 
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Figure 7. Layout of A-Bay. 
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Figure 8. Dimensions (in mm) of Bell Nozzle Test Configuration. 
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Figure 9. Details of Bell Nozzle Test Configuration. Dimensions in mm. *Includes 

Thickness of Gasket. 


Figure 10. Dimensions (in mm) of Stratford Nozzle Test Configuration. 
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2.3. Test Data Systems 

2.3.1. Pressure Measurements 

The static and total pressures were recorded on a Pressure Systems Inc. (PSI) 
Model 780B pressure scanner. The PSI hardware scanned the pressures faster than 
1 00Hz. The PSI hardware averaged the eight most recent scans of each channel then 
transmitted ‘frames’ of all of the channels at slightly less than one frame per second. The 
pressure measurements were recorded four times a second in the ‘facility data’. The 
facility data contains other pressure and temperature measurements of interest to the 
facility operators. Facility data was acquired on a National Instruments PC based system. 

Five, 10 second slices of the facility data were averaged to create five sets of 
nozzle wall static pressure data for the Bell nozzle. These sets of wall static pressure data 
are presented in Section 3.3. The locations of the measurements are described in Section 
2.1 and Appendix 1. The Bell nozzle had about 40 static pressure measurements. The 
Stratford nozzle had no static pressures measured in or on the nozzle. 

2.3.2. Particle Image Velocimetry System 

Dr. Ramachandran, of BAE SYSTEMS Analytical Solutions Inc®, assisted in 
scoping the initial PIV requirements and identifying a company that could perform the 
PIV measurements. BAE subcontracted the PIV measurement task to Dantec Dynamics 
Inc®. Dr. Hammad was the Dantec Dynamics representative responsible for making the 
PIV measurements. Dr. Hammad was involved in the pre-test discussions and made a 
pre-test consultation visit to MSFC to set the final PIV requirements. Dr. Hammad 
traveled to MSFC to make the PIV measurements. Drs. Ramachandran and Hammad 
both helped manipulate and interpret the PIV data. 

In PIV, the flow’s velocity vectors are derived from the movement of particles in 
the flow. The basic PIV process is shown in Figure 1 1 2 , and can be described as follows. 
The flow is seeded with light reflecting particles. The flow in the target area is then 
illuminated with a pulsed laser light sheet. A camera lens images the target area onto the 
charged coupled device (CCD) array of a digital camera. The CCD captures each light 
pulse as a separate image frame. Pairs of images are recorded with a known time 
difference between them. These images are divided into small sub-sections called 
interrogation areas. The interrogation areas from each image frame are cross-correlated, 
pixel by pixel. The correlation produces a signal peak, identifying the common particle 
displacement, dX, for each interrogation area. An accurate measure of the displacement, 
and thus the velocity, is achieved with sub-pixel interpolation. A velocity vector map 
over the whole target area is obtained by repeating the cross-correlation for each 
interrogation area within the image pair. 

The number of velocity vectors in the vector map is dependent upon the number 
of pixels in the CCD camera and the size of the interrogation area. For these entrainment 
tests a CCD camera with 1 186 by 1600 pixels was used for PIV image acquisition. The 
interrogation areas selected were 64 by 64 pixels with 50% overlap in each direction. 

The overlapping technique is used for optimizing the information available in the image 
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maps. With the setup used in these entrainment tests, each PIV dataset resulted in a 
vector map of 36 by 49 data points distributed at regular spacing throughout the camera’s 
field of view (FOV). The CCD camera’s FOV was approximately 133 by 180mm. The 
camera was approximately lm from the laser sheet. Multiple PIV datasets were recorded 
for each nozzle tested. See Sections 3.4. 1.1 and 3.4.2. 1 for specifics on the PIV datasets 
recorded. 

To obtain different datasets the PIV camera was positioned by a traverse system. 

It consisted of two Unidex 1 1 Aerotech® linear stages, each with 300mm traverse length, 
configured in an orthogonal two dimensional traverse system. 

A Dantec Dynamics Inc. PIV system, “FlowMap” 3 , was used in these tests. It 
consisted of a dual pulsed laser system, light sheet optics, the system hub, a personal 
computer (PC) with the FlowMap® software, and GigaBit Ethernet communications link 
between the PC and the system hub. The laser was a New Wave Research® Solo Nd-Yag 
high-energy double oscillator pulsed laser. This class IV laser provided light pulses with 
a maximum energy of 120mJ at a wavelength of 532nm. The laser illuminated at a rate 
of 15 pairs of pulses per second. A Dantec Dynamics 80X60 series light sheet optics 
setup was used to generate a pulsing light-sheet for image plane illumination. The typical 
light sheet thickness was 2 to 3mm. 

The FlowMap System Hub executed the measurement task in uninterrupted 
sequences, while the personal computer performed on-line analysis of the data at the 
performance pace of the PC. Communication between the FlowMap System Hub and the 
PC was via a TCP/IP protocol. The GigaBit Ethernet connections were established in a 
peer-to-peer configuration. The FlowMap System Hub also maintained control of the 
auxiliary devices such as cameras and lasers, and receiving or sending trigger signals to 
or from other devices associated with the measurements. 

The flexible connectivity of the FlowMap System Hub, whether through a 
network or in a peer-to-peer connection, was of great benefit here since the application 
PC was located in a control room, where all monitoring of the experiment took place. 

On-line feedback during data acquisition was accomplished via a FlowMap 
System Hub LIFO (last-in first-out) capability. The LEFO on-line capability manages 
data buffers to guarantee the acquisition performance of the system in uninterrupted 
sequences, while still being able to display on-line images and on-line derived analysis, 
such as PIV velocity vector maps. During PIV data collection, the FlowMap System Hub 
supplied the PC with images at a pace giving the PC time to analyze the data according to 
the complexity of the task and the computing capability of the PC. The LIFO on-line 
capability also allowed for the display of the newest information during the data run, so 
recording quality can be monitored on-line at the same time as the images were safely 
stored at full system performance. No lab time was wasted, when the images were being 
recorded, the on-line monitoring capability ensured a successful off-line post-processing. 
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An air operated seeder (PivPart40 from PIVTEC®) was used to generate seed particles for 
the entrainment tests. Olive oil was used as the seed medium and a typical particle size 
of Ifim. The operating pressure of the seeder was set to about 3bar. 
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Figure 11. The Basic PIV Process. 
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2.3.3. Hot Film Anemometry System 

BAE SYSTEMS Analytical Solutions Inc, an MSFC support contractor, 
developed the hot film anemometry system with hardware available at BAE and MSFC. 
The HFA measurements were made with a TSI® Model Intelligent Flow Analyzer 100. 
The hot film probes were model TSI 1260 A- 10. The probes were mounted in a Unidex 
11 Aerotech traverse mechanism. During a data acquisition run, the traverse mechanism 
was capable of 300mm movement in two axes. The traverse mechanism was mounted on 
a rail system and moved manually between runs as necessary. 

The Labview® programming language was used to drive the traverse mechanism, 
tell the Intelligent Flow Analyzer when to acquire each data point and then store the HFA 
data. It also was used to process the HFA data against the probes’ calibrations. The 
probes were calibrated in a TSI calibrator model 1 125. 
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2.4. Nozzle Set Point Repeatability and Data Accuracy 

This section discusses the variability of the nozzle inlet set point conditions for 
both nozzles. It then discusses the accuracy of the static pressure measurements made for 
the Bell nozzle. The accuracy of the PIV measurements is then discussed. This PIV 
accuracy is an ‘uncertainty’ on the value of the velocity assigned to each data point based 
on the image processing technique used in PIV. It is not a total uncertainty which should 
include at least the effect of: 

• nozzle inlet condition set point variation, 

• variations in the ambient conditions, 

• positional errors in setting up and translating the camera. 

A total uncertainty is not calculated for any of the data presented in this report. 
Repeatability of the PIV measurements is assessed in Section 3.4.4. 

2.4.1. Nozzle Inlet Conditions 

The set points for nozzle inlet conditions for the Stratford and Bell nozzles were 
338.9° K at 10.3 and 19.39atm, respectively. After the test, the facility measurement of 
nozzle inlet total pressure (P to tai) was found to be leaking. The corresponding nozzle inlet 
static pressure (Pstatic) did not leak. When testing the Stratford and Bell nozzles, the area 
ratio for the location at which these facility measurements were made was 18.77 and 
15.42, respectively. Using isentropic gas relations, P to tai should equal P sta ti C /0.9993 and 
Pstatic/0.999, respectively. These ratios result in a P to tai approximately 0.015atm higher 
than Pstatic- This small difference is well within the ±0.15% accuracy (or ± 0.051atm) of 
the 0 to 34atm pressure gages used for these nozzle inlet pressures. Hence, for the test 
data and through out this report, the P sta tic measurement was used for nozzle inlet (i.e., 
total) pressure. 

Stratford Nozzle Inlet Conditions for PIV Acquisition 

The variation of the nozzle inlet conditions during PIV data acquisition for the 
Stratford nozzle is shown in Figure 12 and Figure 13. The average value of the Stratford 
nozzle inlet conditions during PIV acquisition for each data set or ‘window’ is provided 
in Table 1. 
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Table 1. Average Stratford Nozzle Inlet Conditions and Ambient Temperature During 

PIV Data Acquisition 


Window 

Pstatic Avg 

Ttotal Avg 

Tambient Avg 


BSH 

(K) 

(K) 

A 

10.36 

338.30 

291 .88 

E 

10.36 

338.94 

291 .28 

B 

10.35 

339.36 

290.80 

F 

10.36 

338.26 

290.06 

C 

10.36 

338.23 

292.26 

G 

10.37 


293.47 

D 

10.37 

339.69 

293.65 

H 

10.37 


293.79 

1 

10.31 

339.14 

290.50 

J 

10.31 

339.15 

290.62 





B180 

10.36 

339.88 

294.31 

FI 80 

10.35 

337.95 

294.36 

Cl 80 

10.35 

337.07 

294.20 

G180 

10.35 


294.56 

D180 

10.35 

338.08 

294.48 

HI 80 

10.35 

336.76 

294.37 


There was very little variation of nozzle inlet pressure. The I and J windows had 
slightly lower nozzle inlet pressure. The nozzle inlet temperature was usually close to 
339°K, but was as high as 342°K and as low as 334°K. The nozzle inlet temperature was 
measured with a type K thermocouple with an accuracy of ±2.2°K. 
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Bell Nozzle Inlet Conditions for PIV Acquisition 

The variation of the nozzle inlet conditions during PIV data acquisition for the 
Bell nozzle is shown in Figure 14 and Figure 15. The average value of the Bell nozzle 
inlet conditions during PIV data acquisition in each window is provided in Table 2. 


Table 2. Average Bell Nozzle Inlet Conditions and Ambient Temperature During PIV 

Data Acquisition. 


Window 

Pstatic Avg 

Ttotal Avg 

Tambient Avg 


(Atm) 

(K) 

(K) 

A 

19.67 

340.0 

291.6 

E 

19.70 

338.5 

291.8 

B 

19.68 

338.0 

292.5 

F 

19.65 

338.8 

292.5 

C 

19.65 

335.8 

294.0 

G 

19.64 

338.3 

294.1 

D 

19.62 

337.9 

294.6 

H 

19.65 

338.5 

294.7 

K 

19.55 

338.9 

294.7 





B180 

19.66 

337.1 

296.0 

FI 80 

19.68 

339.0 

296.0 

G180 

19.61 

338.4 

296.1 

Cl 80 

19.40 

338.3 

296.1 

D180 

19.19 

337.3 

296.1 

HI 80 

18.92 

337.2 

296.1 


The saw tooth variation in nozzle inlet pressure was relatively small compared to 
its magnitude. During acquisition of the PIV measurements in windows A through G180 
the difference between the maximum and minimum nozzle inlet pressure was about 
0.3atm, or 1.5% of intended set point. This variation was due to the nozzle set point 
being near the NTF maximum supply pressure. As NTF supply pressure dropped, the 
facility valve settings were adjusted to hold temperature and pressure set points close to 
the intended values. 

For windows B 180, FI 80 and G180, the set point was well maintained. Windows 
Cl 80, D180 and HI 80 were obtained as the nozzle inlet pressure began to fall off. No 
effect was seen on the flow entrained in the PIV data . The nozzle inlet temperature only 
varied slightly. The nozzle inlet temperature was measured with a type K thermocouple 
with an accuracy of ±2.2°K. 
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UFA Nozzle Set Points 

The intended nozzle set points for the HFA data acquisition runs were the same as 
the PIV data acquisition runs. The nozzle inlet conditions maintained during the HFA 
runs were very similar to those shown above for the PIV runs. 

2.4.2. Static Pressure Data 

All nozzle test article static pressures were acquired on a Pressure Systems, Inc. 
(PSI) Model 780B pressure scanner. The specifics of the data system are discussed in 
Section 2.3.1. All PSI channels were referenced to atmosphere, except nozzle inlet P to tai- 
Nozzle inlet P s tatic was acquired on a 34atm module and P to tai was acquired on the latm- 
differential module but referenced to nozzle inlet P sta tic instead of atmosphere. 

All pressure measurements made on the PSI hardware had a manufacturer quoted 
uncertainty of 0.15% of the full scale output of the modules they were measured on. PSI 
module noise levels are periodically verified to be within the manufacturer quoted values. 
The modules used for each static pressure measurement are given below. 

• Nozzle Base 1 A through Nozzle Base 2C and Nozzle Inlet P to tai were acquired on 
latm-differential modules. 

• PS1-SR-6 through PSI -SR- 18 & PS1-SR-24 through PS1-SR-36 were acquired 
on 3atm-differential modules. 

• PS2-SR-01 through PS2-SR-05 and PS2-SR-19 through PS2-SR-23 were 
acquired on 6.8atm-differential modules. 

• PS3-SR-37 through PS3-SR-39 and P sta tic were acquired on 34atm-differential 
modules. 
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2.4.3. Particle Image Velocimetry Accuracy 


The method used to estimate the accuracy of the PIV velocity measurements is 
discussed below. The repeatability of the PIV measurements is addressed in Section 


3.4.4. 


To estimate the accuracy of the PIV velocity measurements many parameters 
have to be considered. For the measurements taken in this test, the dominant source of 
error was non-systematic and was mainly due to the uncertainty in the determination of 
the average particle displacement in the interrogation area. A validated method of 
evaluating PIV data uncertainty from laboratory or ‘ideal’ test conditions is provided in 
Reference 2. Paraphrased, in very brief form, it is as follows. 

If a(Axj) is the position uncertainty associated with each individual particle pair, 
then the accuracy of one particle’s displacement estimate originating from an 
interrogation area containing N particles is 

(T(M i ) = jfyl<T 2 (Ax i ) + <r 2 (to2) + "-= -^=T E ^- L 

In practical situations the camera’s F-number is often reduced to a minimum to 
have sufficient light to detect the particles. As a result of the particle diameter, the 
wavelength of light, the camera’s scale factor and F-number, the particle images are less 
than 1 pixel in diameter. If tracking one particle less than a pixel in size, it can be said 
the accuracy is ±0.5 pixel, or approximately that o(Ax) = 0.33 pixels. As a rule of thumb, 
10 particles in the interrogation area are required to have a good signal. With 10 
particles, N=10 and equation 1 results in a position uncertainty of about 0.1 pixels. 

An interrogation area’s average particle displacement in the X and Y directions is 
measured in pixels, Xpix and Ypix, respectively. The magnitude of the particles 
displacement is 


displacement = ^{Xpix 2 + Ypix 2 ) Eq. 2. 

The effect on PIV accuracy of the 0.1 pixels positional uncertainty on each 
interrogation area’s average particle displacement is calculated as follows. In PIV 
measurement terminology, this ‘accuracy’ is called ‘uncertainty’. 

Accuracy = Uncertainty = 0- y n v Eq. 3. 

/ ^(Xpix 2 + Ypix 2 ) 

Or, in terms of a percentage uncertainty, 
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%Uncertainty =100* 


0 . 1 / 


<J(Xpix 2 +Ypix 2 ) 


Eq. 4. 


In PIV the U and V velocity components of the flow are directly scalable from the 
Xpix and Ypix. Therefore, equation 4 can also be applied to quote a percent uncertainty 
for each velocity measurement (that is, at each interrogation area) made by the PIV 
system, 


For these entrainment tests, when processing the data, an interrogation area of 64 
x 64 pixels was selected to maximize the number of seeding particles in an interrogation 
area. Assuming an average of 10 particles per interrogation area would lead to an 
average displacement uncertainty for each interrogation area of about 0.1 pixels. 

These entrainment tests did not present the laboratory or ‘ideal’ test conditions for 
PIV measurements that this 0.1 pixels positional uncertainty was based on. To 
compensate for this, Dr. Hammad, the Dantec Dynamics PIV expert, suggested using a 
more conservative 0.2 pixels average displacement uncertainty. Therefore, for these 
entrainment tests, an uncertainty (for each interrogation area) was calculated to gage the 
PIV data accuracy as follows: 


f 


\ 


%Uncertainty = 100 * 



Eq. 5 


The use of this calculated PIV uncertainty is discussed in Sections 3.4. 1.2.2 and 
3.4.2.2.2 of this report. 
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2.4.4. Hot Film Anemometry Uncertainty 

Two HFA probes were used and both were calibrated for two velocity ranges: a 
‘low speed’ calibration and a ‘high speed’ calibration. For the low speed calibration a 
fourth-order polynomial was utilized to fit the data. Typical accuracy is of the order of 
±2% for speeds more than lm/s. The accuracy decreases below lm/s. At 0.1 m/s the 
accuracy is ±8%. 

For the high speed calibration the data was linearized using the King’s law 
approach as implemented by Collis and Williams. 6 This approach has a maximum error 
of -15% at the higher velocities with an average error around 8%. 7 

The first probe was used for the Stratford nozzle 180° and Bell nozzle 0 and 180° 
planes. Its low speed calibration range was 0.5 to 2.8m/s. Its high speed calibration 
range was 18.0 to 88.0m/s. Due to the limitations of the calibration hardware, there was 
no explicit calibration data for the 2.8 to 18.0m/s range velocities. However, the high 
speed calibration data could be applied all the way down to 0.5m/s and produce velocities 
less than 10% different from those calculated with the low speed calibration. Therefore, 
the low speed calibration was applied for velocities up to 3.0m/s and the high speed 
calibration was applied for velocities above 3.0m/s. 

The second probe was used for the Stratford nozzle 0 and 90° and the Bell nozzle 
90° planes. Post test data analysis determined its low speed calibration was suspect. 
Therefore, its high speed calibration was used for the full velocity range. Based on the 
comparison of the first probe’s low and high speed calibrations the error introduced in the 
low speed range should be less than ±10%. 

Another important factor in the measurement accuracy is the angular response of 
the hot film. The single wire probes used here are most accurate when the probe wire is 
normal to the velocity. If the velocity is at a 45° angle to the wire, the tangential 
component introduces an error, causing the reading to be higher. For regions of highly 
angled flow relative to the wire, the error of the component normal to the wire is 
estimated to be about 15% (high). 

Although no rigorous stack up of the uncertainties was done the uncertainty for 
the HFA data, it could be as much as -10 to +15% for the HFA data. However, this is 
believed to be a conservative uncertainty. 
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3. TEST RESULTS 

The experimental data acquired is discussed in this section. The data acquired 
consists of: schlieren images of the plume structure, laser and fog flow visualization for 
the Bell nozzle, wall static pressures for the Bell nozzle, particle image velocimetry for 
the Stratford and Bell nozzles and hot film anemometry for the Stratford and Bell 
nozzles. 


3.1. Schlieren Images 

Schlieren images (Figure 16) were recorded of the Stratford and Bell nozzle 
plumes prior to modifying the test facility. The test cabin portholes were removed for 
these tests so that the test cabin pressure would remain close to ambient pressure. The 
Stratford nozzle plume is shown for nozzle inlet pressure of 10.24 atmospheres and a 
cabin pressure of 0.977 atmosphere. The Bell nozzle plume is shown for inlet pressure of 
20.06 atmospheres and a cabin pressure of 0.978 atmosphere. The Bell nozzle did not 
quite flow full at this nozzle pressure ratio. The nozzle flow separated from the nozzle 
wall just upstream of the nozzle exit. A Mach disk is visible just downstream of the 
nozzle exit plane. 



Figure 16. Schlieren Images of Stratford and Bell Nozzle Plumes (Flow Right to Left). 
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3.2. Flow Visualization 

Flow visualization is provided for the Bell nozzle. Still images from the flow 
visualization video are presented for three regions: first in the 0 and 1 80° planes - aligned 
with the plume flow; second, normal to the plume flow, just downstream of the nozzle 
exit; third, underneath the plume capture pipe. Two videos of flow visualization, for the 
second and third regions listed above, are included in the electronic media associated 
with this report. 


Aligned With the Plume 

Figure 17 shows the Bell nozzle with the PIV laser illuminating the seed particles 
in the entrained flow. These images, recorded with the facility video camera, are from a 
test while the seeding and laser setup were still being refined. The nozzle inlet pressure 
was approximately lOatm, not at the set point of 19.65atm. The entrainment over the lip 
of the nozzle, into the nozzle, is visible in the 0° plane. 


Normal to the Plume 

Figure 18 shows the hardware orientation used to record the flow visualization in 
a plane normal to the plume. An alternate, lower power laser was used for this flow 
visualization. The laser light sheet generator is visible in Figure 18, pointing south, 
across the plume. The laser sheet was approximately 100mm downstream of the nozzle 
exit. The fog generator was positioned at different times, to put fog both in the base area 
upstream of the nozzle (the test cabin west face) and above the test cabin. The camera 
was looking upstream. 

Figure 19 is a still image from the video showing the nozzle, with fog going into 
the base area and the laser sheet illuminating the fog downstream of the nozzle. The Bell 
nozzle was at its set point of approximately 19.4atm. Figure 20 is nine still images from 
the video which presents about 0.25 seconds of the video. Filaments of entrained flow 
are visible. They were not symmetric, nor were they steady. The flow entrained from the 
base area tended to roll up into large vortices. The vortices were generally attached to the 
west face of the test cabin. The ring in the middle of the images is that fog which had 
mixed into the plume shear layer. 

The video these images were taken from is five minutes long, approximately 1Gb 
in size and is named ‘BellFlowViz.mov’. 


Plume Capture Pipe 

Figure 21 shows the hardware orientation used to record the flow visualization of the 
flow under the plume capture pipe. The flow visualization was done to illustrate the 
vortex in that region. The laser sheet was horizontal, about 100mm below the pipe. 
Figure 22 and Figure 23 are stills from the video. They show the vortex below the pipe 
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inlet. In Figure 22 the overhead lights are still on. The laser’s illumination of the fog and 
vortex is visible. This vortex was attached to the floor, moved around quite a bit, but was 
generally just below the pipe. No other vortex was seen associated with the inlet of the 
plume capture pipe. 

i 

The vortex was large enough at times that it affected the downstream PIV and 
HFA measurements in the 180° plane. Its effect can be seen in the lower right comer of 
the lower image in Figure 17. That region has no fog because the vortex has pushed the 
fog out of the plane of the PIV laser. 

The video these still images were taken from is two minutes long, 444Mb in size 
and is named ‘CapturePipeFlowViz.mov’. 
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Figure 18. Setup for Flow Visualization Normal to Plume. 



Figure 19. Flow Visualization Normal to the Plume, With Overhead Lights On. 
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Figure 21. Setup for Flow Visualization of Flow Under the Plume Capture Pipe. 



Figure 22. Flow Visualization of Flow Under the Plume Capture Pipe, With Overhead 

Lights On. 
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Figure 23. Two Still Images from Flow Visualization of Vortex Under the Plume Capture Pipe. 
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3.3. Nozzle Wall Static Pressure Data 

Figure 24 plots normalized Bell nozzle wall pressures for five nozzle pressure 
ratios (NPR), Piniet/Pambiem- The data plotted in Figure 24 is provided in Table A4. Three 
of the NPRs correspond to a Pi n i et of 1 9.65atm which was the average Pj„iet during PIV 
data acquisition. The Bell nozzle did not flow full at the 19.65atm. The flow separated 
from the nozzle about 30mm upstream of the nozzle exit. 

The pressure ports at X= 146mm were at the nozzle exit plane, on an aft facing 
surface. The pressure recorded by those ports was close to ambient pressure. Note that 
the wall pressures in the nozzle separation were below ambient. These sub-ambient 
pressures induced a region of reverse flow into the nozzle. The flow entrainment near the 
nozzle exit was enhanced by this suction of ambient air into the nozzle. 

The NPR of 19.15 corresponds to the Pj n iet during PIV data acquisition of the 
HI 80 window. The nozzle flow was separated slightly further upstream. The NPR 10.41 
corresponds to Pj n i e t of 10.27atm. This NPR was included for reference only: no 
entrained flowfield data was acquired for the Bell nozzle at this NPR. 



Axial Station Relative to Throat (mm) 

Figure 24. Normalized Bell Nozzle Wall Pressures. 
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3.4. Particle Image Velocimetry Velocity Data 

The PIV data for the Stratford nozzle is discussed first. The PIV data for the Bell 
nozzle is then discussed. The order and format in which the data are presented is very 
similar for both nozzles and is as follows. First, notes about the test data and the relative 
positions of the ‘windows’ are discussed. The velocity field is then discussed in two 
formats; ‘best data without manipulation’ and ‘merged data files’. The first presents the 
best test data for the individual windows from the different delta-T datasets recorded for 
each window. The ‘merged data files’ section discusses the merging of all data in each 
plane into one dataset based on the PIV data uncertainty. It is this merged data that most 
users will find useful. 

The velocity data is then presented in line plots. After discussion of the Stratford 
and Bell nozzle data individually, the two flowfields are compared. 

3.4.1. Stratford Nozzle 

3. 4. 1.1. Data Taken 

Chronological Order with Notes from Run Logs 

All the PIV data could not be taken in one run. Typically two to four windows of 
PIV data could be obtained in each run. The limiting factors were the position of the PIV 
laser sheet generator, the traversing limits of the mechanism that held the PIV camera and 
the NTF’s supply of hot air. The Stratford nozzle PIV data were taken in the following 
order. 

Run 40, Windows A and E 

Run 41, Windows B and F 

Run 42, Window C 

Run 43, Window G 

Run 44, Windows D and H 

Run 45, Windows B 180, FI 80 and Cl 80 

Run 46, Windows G180, D180 and HI 80 

Run 49, Windows I and J. 

For the PIV technique to work well, the test medium must be homogeneously 
seeded. The test medium for this entrainment experiment was the air in (and that air 
flowing into) the A-bay of building 4777. For the first several runs the seeding was 
either insufficient or poorly diffused. This marginal seeding affected windows A, E, B, F 
and C. Several different combinations of seeders and seeding locations were tried to 
correct the poor seeding. By Run 43, window G, a good location was found for the 
seeder. Another factor influencing the seeding quality was the window in the south wall 
(see Section 2.2.2) of the A-bay. This window was covered with a loose fitting tarp 
(starting with window G), significantly reducing the air flow through it. At that point the 
quality of the seeding was much improved and difference in the quality of the PIV data 
after these changes is visible in the contour plots shown in the following section, Section 

3.4. 1.2. 
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For Run 45, the rollup door on the west end of building 4777 was not opened. 
This was an oversight. This does not appear to have affected the data in windows B1 80, 
F180 and C180. 

Delta-T Interval Matrix and File Naming Convention 

The term 'window’ refers to a region in space in which data was collected. This 
region was function of the camera’s field of view (FOV). The term ‘dataset’ refers to a 
set of PIV data recorded at a specific delta-T in a window. Multiple PIV datasets were 
obtained at each window. Table 3 documents the different delta-T intervals used to 
acquire the PIV datasets. Each of these PIV datasets is in the electronic media associated 
with this report. 


Table 3. Matrix of Delta-Ts at Each Window and Names of the Data 

Files. 


Laser PI use Dt 


WINDOW 

2 ms 

5 ms 

20 ms 



200 ms 


500 ms 

A 




Sonic A050.dat 

Sonic A100.dat 

Sonic A200.dat 


Sonic A500jdat 

B 



^1'. 


Sonic B100.dat 




C 


Sonic C005.dat 


Sonic C050.dat 

Sonic C100.dat 

Sonic C200.dat 


Sonic C500.dat 

D 


Sonic D005.dat 


Sonic D050.dat 

Sonic D100.dat 

Sonic D200.dat 


Sonic D500.dat 

E 









F 





Sonic F100.dat 

Sonic F200.dat 


Sonic F500.dat 

G 




Sonic G050.dat 

Sonic G100.dat 

Sonic G20Q.dat 


Sonic G500.dat 

H 

■ 

■1 V 






Sonic H500.dat 


■■ 

■i ■ 







B180 


Sonic B180005.dat 


Sonic B180050.dat 

Sonic Bl80100jdat 

Sonic B180200.dat 


Sonic B18050adat 

C180 





Sonic Cl80100dat 

Sonic Cl80200.dat 


Sonic Cl80500.dat 

D180 

Sonic D180002.dat 

Sonic D180005idat 


Sonic D180050.dat 

Sonic Dl80100dat 

Sonic D180200.dat 


Sonic D180500.dat 

F180 





Sonic F180100.dat 

Sonic F180200.dat 


Sonic F180500.dat 

G180 

■ ■! 

!■ 

1 ■ 






H180 





<0 

■q 

o 

o 

o 

00 

X 

o 

c 

o 

CO 

Sonic H180200.dat 


Sonic HI 8050adat| 




i«aL>K0Hi 


Sonic I100.dat 

Sonic I200.dat 


Sonic ISOQdat 1 










Relative and Absolute Positions of PIV Windows 

A schematic of the windows in which PIV datasets were recorded is shown in 
Figure 25. Eight windows of PIV data were obtained above the nozzle in the 0° plane. 
Two windows were obtained above the plume capture pipe. Six windows of PIV data 
were obtained below the nozzle in the 180° plane. 

The position of the nozzle in the camera’s FOV in the A and B and B 1 80 
windows was used to determine the coordinate transformations required to locate 
windows A and B 180 in the nozzle global coordinate system of the Stratford nozzle. 
(Note that the Stratford and Bell nozzle global coordinate systems are different). The 
data from the rest of the windows in the 0° and 1 80° planes were then translated relative 
to the A and B180 windows based on the commands sent to the traversing mechanism. 

The position of the plume capture pipe in the camera’s FOV in the I window was 
used to determine the coordinate transformations required to locate I in the coordinate 
system of the pipe. Window J’s data was then translated relative to I based on the 
commands sent to the traversing mechanism. 


- 42 - 













FinalReport8.2 


11/15/2005 


The translations for all windows are provided in Table 4. These translations were 
applied to the coordinate system of the datasets for each window. Each window’s 
datasets consist of a matrix of 36 by 49 data points. The origin of the coordinate system 
of each window’s data points is that window’s lower left hand comer. The fifth and sixth 
columns in each dataset file (Table 3) are the X and Y positions (in mm) of each data 
point within that window’s datasets. The translations given in Table 4 were applied to 
these X-Y coordinates to position each window’s datasets in the Stratford nozzle’s global 
coordinate system. 

The coordinate systems within the datasets in the 0° plane are identical. That is, 
the matrix of X-Y locations within the datasets relative to the windows’ origin is the 
same. Likewise all datasets in the 1 80° plane have the same matrix of X-Y points. There 
are slight differences between the 0° and 1 80° planes. See notes in Table 4 for details. 
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3.4. 1.2. Velocity Field Discussion 

The velocity data is presented in two forms. First, the best dataset for each 
window, without modification of the datasets, is presented. Second, a single dataset 
created from the best data available from datasets of all windows, is presented. 

3.4. 1.2.1. Best Data Without Manipulation 

Table 3 lists the different delta-Ts used to acquire data for each PIV window. The 
smaller delta-T’s resolved higher velocities better than the larger delta-T’s. Figure 26 
illustrates this. The 50 microsecond delta-T resolved higher velocities close to the nozzle 
exit and in the plume shear layer. However, the smaller delta-Ts have higher uncertainty 
in the low velocity regions as illustrated in Figure 27. 

Of the data recorded, it was judged that the best representation of the flowfield 
(without modification of the datasets) is captured by the 100 microsecond delta-T dataset 
for windows A to D and B180 to D180 and 200 microsecond delta-T dataset for windows 
E to K and FI 80 to HI 80. Figure 28 presents these datasets. Note in these TecPlot® 
images that the overlapped regions of adjacent windows are plotted, but only the window 
that is in ‘on top’ is visible. 

The datasets in Figure 28 are the best compromise between capturing the high 
velocity regions, acceptable uncertainty levels and smooth velocity gradients between 
adjacent windows. These datasets show the high velocity near the nozzle exit and the 
general acceleration of the entrained flow with axial position. The seeding was less than 
optimum for windows A, B, C, E and F and the resulting PIV data is of reduced quality 
for those windows. Window C was particularly difficult to get good seeding and that 
data has obvious discontinuities with its neighboring windows. The seeding was 
adequate for windows D, G, and H in the 0° plane and for all the 1 80° windows. The PIV 
data in those windows appears to be better quality data. 

The character of the flow is symmetric. The PIV data indicates the 1 80° plane 
had higher peak velocities near the nozzle exit. This was probably caused more by the 
poor seeding in the A and B windows than by real physics. 

The indicated low velocity near the centerline is the plume flow. The velocity 
was not measurable with the delta-Ts presented here. Also, the plume was not seeded as 
the intent of the test was to measure the entrained flow. The apparent low velocity region 
at the outer edge of E to H is discussed in Section 3.4.3. 


- 46 - 



FinalReport8.2 


11/15/2005 


Sonic Odeg Plane 



X (mm) 


Figure 26. Velocities Resolved with Different Delta-Ts. 
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Figure 27. Uncertainty Resulting from Different Delta-Ts. 
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Figure 28. Best Representation of Stratford Nozzle Entrainment Flowfield without 

Modification of the PIV Datasets. 
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3.4.1 .2.2. Merged Data Files 

Velocity Contour Plots 

A more accurate representation of the entrained flow was created by combining 
the data from the multiple delta-T datasets available for each window. The high velocity 
regions from the small delta-T datasets were ‘merged’ into the larger delta-T datasets that 
captured the general entrained flowfield better. Merging of multiple delta-T datasets was 
only required for the inner row of windows (A to D and B 1 80 to D 1 80) where there were 
large gradients in measured velocity. 

Figure 26 and Figure 27 are used to explain the merging of multiple delta-T 
datasets. These two figures show that the 50 microsecond delta-T resolved higher 
velocity near the nozzle exit and in the nozzle shear layer with relatively low 
uncertainties. If one were willing to accept 8% uncertainty (Figure 27), much of the 50 
microsecond delta-T dataset is acceptable near the plume shear layer and nozzle exit. 
These higher velocities were incorporated into the merged data shown in Figure 29. 

The merging was done by assessing the data, point by point. The logic for 
merging was that the data from the smallest delta-T dataset was used if the uncertainty at 
that data point was ‘acceptable’ - in this example - an uncertainty of less than 8%. 
Because the calculation of the PIV data uncertainty is somewhat subjective (as described 
in Section 2.4.3), an uncertainty was chosen that resulted in a fairly continuous velocity 
flowfield. For the Stratford nozzle data, 8% uncertainty was deemed a good compromise 
between accuracy and resolution of the flowfield. Lower uncertainty levels (2 and 5% 
were assessed) resulted in noticeable flowfield discontinuities or flowfields that failed to 
capture the higher velocities near the plume shear layer and nozzle exit. 

The merging of data within each window was done in an Excel spreadsheet. 

These Excel files are provided in the electronic media associated with this report. The 
structure of these files is discussed in Section 3.4.6. If an alternate uncertainty is desired, 
these files can be manipulated to produce merged data with the alternate uncertainty as 
the criterion. 

After the data was merged within each window, the merged data were then joined 
or ‘woven’ together (by sorting on the radial and axial stations in the spreadsheet). Thus, 
creating a single data field that encompassed all of windows A to F (Figure 29) and a 
second that encompassed windows B 1 80 to H 1 80 (Figure 31). These woven data fields 
allow the overlapped regions to be seen when plotting contours in TecPlot®. They also 
enable extracting line plots through the entire flowfield. (See Section 3.4.6 for discussion 
of the Excel files). 

To clarify the terminology used here: ‘merged’ means that data from multiple 
delta-T datasets were evaluated point by point to determine the best value to use at each 
data point in a window. ‘Woven’ means that the data points from multiple windows were 
reordered into one contiguous data field for plotting. 
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Figure 29 and Figure 3 1 show that the resulting merged and woven PIV data 
resulted in fairly continuous flowfields. In these velocity contours, the saw tooth patterns 
indicate velocity discontinuities in the overlapped regions of adjacent windows. The 
relatively poor quality of the measurements in the C window is evident by the saw tooth 
patterns surrounding it. The 1 80° plane had much better agreement window to window 
and in general appears to be a better set of data than the 0° plane. 

Figure 30 and Figure 32 show the corresponding uncertainties. In the uncertainty 
plots, the strong gradients between the inner and outer rows of windows are due to the 
different delta-Ts and are not a result of velocity discontinuities. Note the uncertainties 
near the centerline are acceptable for windows D, B 1 80, C 1 80 and D 1 80. This is because 
those windows had delta-Ts of 2 and 5 microseconds and sufficient seed particles had 
diffused into the plume shear layer. Also note the band of greater than 8% uncertainty (in 
Figure 30) along the top of the 0° plane data. This phenomenon is related to the apparent 
low velocities in this same area (see Figure 29) and is discussed in Section 3.4.3. 

Merged and woven data for the 0° and 1 80° planes are shown together in Figure 
33. Again, the general character of the 0° and 180° planes agree. The 180° plane appears 
more continuous than the 0° plane and indicates higher velocities near the nozzle exit. 

The contours in Figure 33 show there were two regions where the entrained flow 
was accelerated. These regions of acceleration indicate the location of the flow 
phenomena that induced the entrained flow. There was a distinct acceleration of the 
entrained flow near the nozzle exit, and a second more gradual acceleration, with axial 
station, of the entire entrained flowfield. 

Figure 34 shows the velocity vectors in the flowfield near the nozzle exit. The 
vectors are not scaled by their magnitude. 
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Data merged in A to D for 8% Uncertainty . 
Data woven together for A to H . 
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Figure 29. Merged and Woven PIV Velocity for 0° Plane for the Stratford Nozzle. 
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Figure 30. Uncertainty for the Merged and Woven PIV Data for 0° Plane for the 

Stratford Nozzle. 
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Data woven together for B1 80 to H 1 80. 


Figure 3 1 . Merged and Woven PIV Velocity for 1 80° Plane for the Stratford Nozzle. 
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Figure 32. Uncertainties for Merged and Woven PIV Data for 180° Plane for the 

Stratford Nozzle. 
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Figure 33. Merged and Woven PIV Velocity for Stratford Nozzle. 
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Figure 34. Velocity Vectors from the Merged and Woven PIV Data for the Stratford 

Nozzle. 
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Velocity Line Plots 

Line plots of velocities at constant radial and axial stations are now discussed to 
illustrate some of the flowfield characteristics. The radial and axial stations for these line 
plots are shown in Figure 35. These plots were created in the Excel file that generated 
the woven data. The data plotted here can be extracted from that Excel file. The Excel 
file can be used to create additional line plots at alternate radial or axial stations as well. 

Figure 36 and Figure 37 show the velocities at five constant radial stations from 
the merged and woven data in the 0° and 1 80° planes. The saw tooth patterns in the lines 
are the overlapping regions of adjacent windows. In Figure 36 the C window data is 
inconsistent with its neighbors and should be disregarded. The Cl 80 window data in 
Figure 37 indicates the velocity plots should be smooth across windows B, C and D. The 
zero velocities in window A of Figure 36 is where those radial stations intersect the 
nozzle hardware. 

These velocity line plots show, as did the velocity contour plots, that there were 
two regions where the entrained flow was accelerated. The first region was the strong 
acceleration close to the nozzle exit. The flow phenomenon inducing this acceleration 
was the viscous pumping of the plume shear layer. The second region of acceleration 
was the more general acceleration of the entire entrained flowfield with axial station. 

This second acceleration started at about X= 125mm. The primary flow phenomena 
inducing this second acceleration was the sink created by the interaction of the plume and 
the plume capture pipe. The cause of the plume capture pipe sink is discussed in Section 
3.4. 1.2. 3. 

The peak magnitudes shown in these line plots are approximately, 1 1 .5 and 13m/s 
for the 0 and 1 80° planes. The U-velocity component peaked about 25mm downstream 
of the nozzle exit (X=0mm). The V-velocity component peaked about 50mm 
downstream of the nozzle exit. The nozzle outer radius at the exit was 20.3mm. The 
acceleration near the nozzle exit is evident in the profiles out to about 134mm 
(Y= 134mm profile not shown here), or about 1 14mm above the nozzle lip. The effect of 
the nozzle exit acceleration extended about 125mm downstream of the nozzle exit. 
Downstream of X~125mm the entrained flow’s velocity increased due (mainly) to the 
plume capture pipe sink’s effect. 

The erratic and decreasing values of the Y=44mm profiles starting at X= 100mm 
indicates bad data. In this region, the PIV measurements were affected by the shear 
layer. 

Figure 38 compares three profiles for the 0 and 180° planes. They compare quite 
well considering the less than satisfactory seeding quality in the first windows of the 0° 
plane measurements. 

Figure 39 and Figure 40 present the 0° plane velocities at constant axial stations. 
Figure 41 presents the 180° plane velocities at constant axial stations. In Figure 39, the 
peak velocities near the nozzle exit are illustrated. Note that the X=- 100.4, -39.1 and - 
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1 .2mm profiles intersect the nozzle hardware at their lower ends. Hence, the velocities 
go to zero for these profiles. The X=0.6 and 20.5mm profiles go toward zero due to the 
absence of seed particles in the plume. The X=-1.2 and 0.6mm profiles indicate there is 
not a significant difference within the overlapped region of windows A and B. As one 
would expect, the peak velocities increased closer to the nozzle exit. 

The X=200 and ~300mm velocity profiles in Figure 40 and Figure 41 illustrate 
the acceleration of the bulk of the entrained flow with increasing axial station. Again, 
this acceleration was primarily due to the plume capture pipe sink. In Figure 40, the poor 
seeding in the C window was avoided by plotting the X=1 98.9mm profile from window 
D data. 


Note that the velocities decreased, or trended toward zero, near the centerline - 
that part of the profile closest to or in the plume shear layer. One would expect the flow 
to be accelerated by the high velocity shear layer. However, the plume was not seeded. 
Therefore, any PIV measurements of the shear layer would have required some of the 
seed material to have diffused from the ambient air into the plume shear layer. In 
addition, smaller delta-Ts would be required to adequately resolve this high speed shear 
layer flow. Therefore, the velocities trending to zero in this region should be disregarded. 

For the profiles between X=40 and 300mm the peak velocities occurred further 
from the nozzle centerline with increasing axial station. This was due to the increasing 
diameter of the plume shear layer with axial station. 

The 1 80° plane velocities have the same general trends as the 0° plane. Figure 42 
compares some of the axial profiles for both planes. 
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Figure 35. Axial and Radial Stations of Subsequent Velocity Profiles for the Stratford 

Nozzle. 
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Figure 36. PIV Measured Velocities at Five Radial Stations in the 0° Plane for the 

Stratford Nozzle. 
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180 Degree Plane, 8% Uncertainty 
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Figure 37. PIV Measured Velocities at Five Radial Stations in the 180° Plane for the 

Stratford Nozzle. 
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0 and 180 Degree Planes, 8% Uncertainty 
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Figure 38. Comparing PIV Measured Velocities for the 0 and 180° Planes at Three 
Radial Stations for the Stratford Nozzle. 
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0 Degree Plane, 8% Uncertainty 
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Figure 39. PIV Measured Velocities at Five Axial Stations in the 0° Plane for the 

Stratford Nozzle. 
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0 Degree Plane, 8% Uncertainty 
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Figure 40. PIV Measured Velocities at Five Additional Axial Stations in the 0° Plane for 

the Stratford Nozzle. 
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180 Degree Plane, 8% Uncertainty 
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Figure 41 . PIV Measured Velocities at Six Axial Stations in the 1 80° Plane for the 

Stratford Nozzle. 
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0 and 180 Degree Planes, 8% Uncertainty 
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Figure 42. Comparing PIV Measured Velocities for the 0 and 180° Planes at Axial 

Stations for the Stratford Nozzle. 
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3. 4. 1.2. 3. Windows I and J - The Plume Capture Pipe Sink 

For the Stratford nozzle, two additional windows of PIV data were acquired just 
above the plume capture pipe. Windows I and J are shown in Figure 43 through Figure 
45. The merged data for the I window was created from datasets with significantly 
different delta-Ts. This explains in the unusual velocity contours above and in front of 
the pipe inlet. The PIV data resolved maximum velocities of 25m/s above the pipe and 
61 m/s below the upper lip of the pipe. 

These velocities are evidence of sub ambient pressure, or a ‘sink’, which existed 
inside the plume capture pipe. This sink was created by the ejector effect of the plume 
attaching to the inner wall of the plume capture pipe. A viscous, high speed flow, in 
close proximity to a wall, creates a low pressure region. This sub ambient pressure acted 
as a pump that sucked ambient air in. The velocity contours and vectors in Figure 45 
show flow heading inward from all locations - a classic example of a sink. This sink 
existed for the full circumference of the plume capture pipe. 

The Excel files of the PIV data for I and J windows are included in the electronic 
media associated with this report. The structure of the Excel files is similar to those for 
the rest of the PIV windows. These Excel files are described in Section 3.4.6. 


- 65 - 



Y, Relative to Nozzle Centerline (mm) 


FinalReport8.2 


11/15/2005 




Data merged in A to D for 8% uncertainty. 
D ata woven together for All. 


Vel Mag (m/s) 

14 
12 
10 
8 
6 
4 
2 
0 


Data merged in B180 to D180 
for 8% uncertainty. 

Data woven together for B1 80 to HI 80. 


' '' I '' '' I '' '' I '' '' I ' 

400 500 600 700 


X, Relative to Nozzle Exit (mm) 


Figure 43. PIV Data Near the Plume Capture Pipe. 
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Figure 44. Close up View of PIV Data Near the Plume Capture Pipe. 
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Figure 45. Velocity Vectors from the PIV Data Near the Plume Capture Pipe. 
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3.4.2. Bell Nozzle 

3.4.2. 1. Data Taken 

Chronological Order with Notes from the Run Logs 

All the PIV data could not be taken in one run. Typically, two to four windows of 
PIV data could be obtained in each run. The limiting factors were the position of the PIV 
laser sheet generator, the traversing limits of the mechanism that held the PIV camera and 
the NTF’s supply of hot air. The Bell nozzle PIV data were taken in the following order. 

Run 50, Windows A and E 
Run 5 1 , Windows B and F 
Run 52, Windows C and G 
Run 53, Windows D, H and K 

Run 55, Windows B 180, FI 80, G180, Cl 80, D180 and HI 80. 

Bell Nozzle inlet conditions during PIV data acquisition are shown in Figure 14 
and Figure 15. Run 55 was the last possible run of the test series. The intent was to get 
all data possible even if the nozzle set point could not be maintained (see Figure 15). The 
nozzle set point was maintained for windows B 1 80, F 1 80 and G 1 80. It began to drift 
down from ~19.65atm to ~19.4atm during acquisition of window Cl 80. The set point 
continued to drift down during acquisition of windows D180 (~19.2atm) and HI 80 
(~18.9atm). 


Delta-T Matrix and File Naming Convention. 

The term ‘window’ refers to a region in space in which data was collected. This 
region was function of the camera’s field of view. The term ‘dataset’ refers to a set of 
PIV data recorded at a specific delta-T in a window. Multiple PIV datasets were obtained 
at each window. Table 5 documents the different delta-T intervals used to acquire the 
PIV datasets. Each of these PIV datasets is in the electronic media associated with this 
report. 


Table 5. Matrix of Delta-Ts at Each Window and Names of the Data Files. 
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Relative and Absolute Positions of PIV Windows 

A schematic of the windows in which PIV datasets were recorded is shown in 
Figure 46. Nine windows were obtained above the nozzle in the 0° plane. Six windows 
of PIV data were obtained below the nozzle in the 1 80° plane. 

The position of the nozzle in the camera’s FOV in the A, B and B1 80 windows 
was used to determine the coordinate transformations required to locate windows A and 
B 180 in the global coordinate system of the Bell nozzle. (Note that the Stratford and Bell 
nozzle global coordinate systems are different). The data from the rest of the windows 
were then translated relative to the A and B 1 80 windows based on the commands sent to 
the traversing mechanism. The translations for each window are provided in Table 6. 
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Table 6. Window Coordinate Transformations for Bell Nozzle PIV Data. 


General Notes 

Global Origin is at nozzle exit on nozzle centerline 
All dimensions are in millimeters 
The nozzle outer lip radii is 81 .28 

Notes for the Odegree Plane 

Windows A through H and K are 133.072 mm wide by 179.8 mm tall 
The origin of A was 1 19.922 to left of nozzle exit plane 


The bottom of A through D was 22.87mm below the nozzle lip. 

E through H were 155 mm above A through D. K was 250 mm above D 


Wndow Name 


Translations 


Total 

Translation 


Global 

Coordinates 


Wndow Name 

K 


X 

Y 

Translations 

-119.922 +110 
+110 +90 

-22.87 +81.28 
+250 

Total 

Translation 

Global 

Coordinates 

190.078 
190.078 to 
323.16 

308.410 
308.41 to 
488.21 




Wndow Name 

B180 

C180 


X 

Y 

X 

Y 

Translations 

-14.700 

-160.3-81.28 

-14.7 +110 

-160.3 -81.28 

Total 

Translation 

-14.700 

-241.580 

95.300 

-241.580 

Global 

Coordinates 

-14.700 to 
117.7 

-241 .58 to - 
62.98 

95.3 to 227.7 

-241 .58 to - 
62.98 


185.300 -241.580 

-241.58 to - 

185.3 to 317.7 62.98 


Wndow Name 


Translations 

Total 

Translation 

Global 

Coordinates 



-160.3-8 

-14.700 

155 

-14.700 

-396.5 

-14.700 to 

-217.98 

117.7 

396.5 



-160.3-81.28- 
14.7+110 + 90 155 


Notes for 180degree Plane 
Wndows B180 through HI 80 are 132.187 mm wide by 178.6 mm tall 
The origin of B180 was 14.7 mm to left of nozzle exit plane 

C180 & G180 were 1 10 mm to right of B180 & F180. D180, H180 were 90 mm to right of C180 & G180. 

The bottom of B 180 through D 180's was 160.3 mm below the nozzle lower lip. 

FI 80 through HI 80 were 155 mm below B 180 through D180. 
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These translations were applied to the coordinate system of the datasets for each 
window. Each window’s datasets consist of a matrix of 36 by 49 data points. The origin 
of the coordinate system of each window’s data points is that window’s lower left hand 
comer. The fifth and sixth columns in each dataset file (Table 5) are the X and Y 
positions (in mm) of each data point within that window’s datasets. The translations 
given in Table 6 were applied to these X-Y coordinates to position each window’s 
datasets in the Bell nozzle’s global coordinate system. 


The internal coordinate systems of the datasets in the 0° plane are identical. That 
is, the matrix of X-Y locations in the datasets relative to the windows’ origin is the same. 
Likewise all datasets in the 180° plane have the same matrix of X-Y points. There are 
slight differences between the 0° and 180° planes. See notes in Table 6 for details. 




■ VT 
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3.4. 2.2. Velocity Field Discussion 

The velocity data is presented in two forms. First, the best dataset for each 
window, without modification of the datasets, is presented. Second, a single dataset 
created from the best data available from all datasets of all windows, is presented. 

3 .4.2.2. 1 . Best Data Without Manipulation 

Table 5 lists the different delta-Ts used to acquire data for each PIV window. The 
smaller delta-Ts resolved higher velocities better than the larger delta-Ts. Figure 47 
illustrates this. The 5 and 50 microsecond delta-Ts resolved the higher velocities close to 
the nozzle exit. However, these small delta-Ts have higher uncertainty in the lower 
velocity regions (Figure 48). 

Of the data recorded, it was judged that the best representation of the flowfield 
(without modification of the datasets) is captured by the 200 microsecond delta-T 
datasets for windows A to D and B 180 to HI 80 and 500 microsecond delta-T datasets for 
windows E to K. Figure 49 presents these datasets. Note in these TecPlot images that 
the overlapped regions of adjacent windows are plotted, but only the window that is ‘on 
top’ is visible. 

The datasets in Figure 49 are the best compromise between capturing the high 
velocity regions, acceptable uncertainty levels and smooth velocity gradients between 
adjacent windows. These datasets show the high velocity near the nozzle exit and the 
general acceleration of the entrained flow with axial position. The character of the flow 
is symmetric, but the 1 80° plane has lower peak velocities. The indicated low velocity 
near the centerline is the plume flow. The plume velocity was not measurable with the 
delta-Ts presented here. Also, the plume was not seeded as the objective of this test was 
to measure the entrained flow. The apparent low velocities at the outer edge of windows 
E to H and F180 toH180 are discussed in Section 3.4.3. 


- 72 - 



FinalReport8.2 


11/15/2005 


Bell Odeg Plane 



X (mm) 

Figure 47. Velocities Resolved with Different Delta-Ts. 
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Figure 48. Uncertainty Resulting from Different Delta-Ts. 
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Figure 49. Best Representation of Bell Nozzle Entrainment Flowfield without 
Modification of the PIV Datasets. 
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3.4. 2.2.2. Merged Data Files 

Velocity Contour Plots 

A more accurate representation of the entrained flow was created by combining 
the data from the multiple delta-T datasets available for each window. The high velocity 
regions from the small delta-T datasets were ‘merged’ into the larger delta-T datasets that 
captured the general entrained flowfield better. Merging of multiple delta-T datasets was 
only required for the inner row of windows (A to D and B1 80 to D1 80) where there were 
large gradients in measured velocity. 

Figure 47 and Figure 48 are used to explain the merging of multiple delta-T 
datasets. Note in the 5 and 50 microsecond datasets that higher velocities were resolved 
near the nozzle exit. The 50 microsecond delta-T dataset had relatively low uncertainty 
in this region. If one were willing to accept 8% uncertainty, much of the 50 microsecond 
delta-T dataset is acceptable near the plume shear layer and nozzle exit. These higher 
velocities were captured in the merged data shown in Figure 50. 

The merging was done by assessing the data, point by point. The logic for 
merging was that the data from the smallest delta-T dataset was used if the uncertainty at 
that data point was ‘acceptable’ - in this example - an uncertainty of less than 8%. 
Because the calculation of PIV data uncertainty is somewhat subjective (as described in 
Section 2.4.3) an uncertainty was chosen that resulted in a fairly continuous velocity 
flowfield. For the Bell nozzle data, 8% uncertainty was deemed a good compromise 
between accuracy and resolution of the flowfield. Lower uncertainty levels (2 and 5% 
were assessed) resulted in noticeable flowfield discontinuities or flowfields that failed to 
capture the higher velocities near the plume shear layer and nozzle exit. 

The merging of data within each window was done in an Excel spreadsheet. 

These Excel files are provided in the electronic media associated with this report. The 
structure of these files is discussed in Section 3.4.6. If an alternate uncertainty is desired, 
these files can be manipulated to produce merged data with the alternate uncertainty as 
the criterion. 

After the data was merged within each window, the merged data were then joined 
or ‘woven’ together (by sorting on the radial and axial stations in the spreadsheet). Thus, 
creating a single data field that encompassed all of windows A to F (Figure 50) and a 
second that encompassed windows B 180 to HI 80 (Figure 52). These woven data fields 
allow the overlapped regions to be seen when plotting contours in TecPlot. They also 
enable extracting line plots through the entire flowfield. (See Section 3.4.6 for discussion 
of the Excel files). 

To clarify the terminology used here: ‘merged’ means that data from multiple 
delta-T datasets were evaluated point by point to determine the best value to use at each 
data point in a window. ‘Woven’ means that the data points from multiple windows were 
reordered into one contiguous data field for plotting. 
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Figure 50 and Figure 52 show that the resulting merged and woven PIV data 
resulted in fairly continuous flowfields. The saw tooth patterns indicate velocity 
discontinuities in the overlapped regions of adjacent windows. Figure 5 1 and Figure 53 
are the corresponding uncertainties. For the uncertainties, the strong gradients between 
the inner and outer rows of windows are due to the different delta-Ts, not flow 
discontinuities. 

Merged and woven for the 0° and 180° planes are shown together in Figure 54. 
Again, the general character of the 0° and 1 80° planes agree, but the peak magnitudes are 
lower in the 1 80° plane. 

The contours in Figure 54 show there were two regions where the entrained flow 
was accelerated. These regions of acceleration indicate the location of the flow 
phenomena that induced the entrained flow. There was a distinct acceleration of the 
entrained flow near the nozzle exit, and a second more gradual acceleration, with axial 
station, of the entire entrained flowfield. 

Figure 55 shows the velocity vectors near the nozzle exit. The vectors are not 
scaled by their magnitude. The effect of the reverse flow into the separated region of the 
Bell nozzle is evident in the vectors just downstream of the nozzle exit. 
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Figure 50. Merged and Woven PIV Velocity for 0° Plane for the Bell Nozzle. 
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Figure 5 1 . Uncertainty for the Merged and Woven PIV Data for 0° Plane for the Bell 

Nozzle. 
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Figure 52. Merged and Woven PIV Velocity for 180° Plane for the Bell Nozzle. 
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Figure 53. Uncertainty for the Merged and Woven PIV Data for 180° Plane for the Bell 

Nozzle. 
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Figure 54. Merged and Woven PIV Velocity for Bell Nozzle. 
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Figure 55. Velocity Vectors from the Merged and Woven PIV Data for Bell Nozzle. 
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Velocity Line Plots 

Line plots of velocities at constant radial and axial stations are now discussed to 
illustrate some of the flowfield characteristics. The radial and axial stations for these line 
plots are shown in Figure 56. These plots were created in the Excel file that generated 
the woven data fields. The data plotted here can be extracted from that Excel file. The 
Excel file can be used to create additional line plots at alternate radial or axial stations as 
well. 


Figure 57 shows the velocities at five constant radial stations from the merged and 
woven data field in the 0° plane. The saw tooth patterns in the lines are the overlapping 
regions of adjacent windows. 

These velocity line plots show, as did the velocity contour plots, there were two 
regions where the entrained flow was accelerated. The first region was the strong 
acceleration close to the nozzle exit. For the Bell nozzle this acceleration was due to two 
phenomena, the nozzle flow separation and the viscous pumping of the plume shear 
layer. The second region of acceleration was the more general acceleration, with axial 
station, of the entire entrained flowfield. This second acceleration started at about 
X= 125mm. The primary flow phenomena inducing the second, more general, 
acceleration was the sink created by the interaction of the plume and the plume capture 
pipe. The cause of the plume capture pipe sink is discussed in Section 3.4. 1.2. 3. 

The acceleration near the nozzle exit is evident in all five profiles shown but most 
dramatically in the Y=94mm profile. The U-velocity component peaked just upstream, 
X~-10mm, of the nozzle exit (X=0mm) and the V- velocity component peaked just 
downstream of it, X~+10mm. The rapid decrease of the U-velocity component was a 
result of the reverse flow into the bell nozzle. The nozzle outer radius at its exit was 
81.28mm. Although not shown in these line plots the effect of the nozzle exit was seen 
out to a radius of Y=~227mm in the 0° plane (and Y=~213mm in the 180° plane). The 
effect of the nozzle exit acceleration extended about 125mm downstream of the nozzle 
exit. 


Downstream of X~125mm, the entrained flow’s velocity increased due to 
(mainly) the plume capture pipe sink’s effect. The slight erratic shape of the Y=94mm 
profile in the C window was caused by the difficulty in getting good PIV seeding for that 
window. 

The same five radial stations are presented for the 180° plane in Figure 58. These 
profiles are not as smooth but exhibit the same general trends as the 0° plane. Figure 59 
compares profiles for the 0 and 1 80° planes. Other than the velocity magnitudes close to 
the nozzle exit they agree well. 

Figure 60 presents the 0° plane velocity at constant axial stations from the left 
hand side of windows A and E to just downstream of the nozzle exit. The peak velocities 
near the nozzle exit are illustrated. Note that the X=-40.8 and -1 .2mm profiles intersect 
the nozzle hardware at their lower ends. Hence, the velocities go to zero for these 
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profiles. The X=-l .2 and 0.9mm profiles agree well and indicate there is not a significant 
difference within the overlapped region of windows A and B. They both show zero 
velocity at their lower ends because of the proximity of the nozzle hardware. Note that 
the 1 9mm profile contains negative U-velocity component. This resulted from the 
reverse flow induced by the separated nozzle flow. 

Figure 61 plots 0° plane velocities at constant axial stations between X=40 and 
300mm. The U-velocity components at X=40.4 and 62.0mm show the affect of the 
reverse flow near the nozzle exit. They both reached a maximum at Y~125mm then 
decreased to a local minimum at Y~75mm. Downstream of X=125mm, the bulk of the 
entrained flow began to accelerate. This acceleration was in part due to the shear layer 
but mainly a result of the plume capture pipe sink. 

For the velocity measurements between X=0 and ~ 100mm the uncertainties (in 
Figure 51) indicate these velocities were valid. Sufficient seeding mixed into the plume 
shear layer in the nozzle separation induced recirculation. However, downstream of 
X= 100mm, the velocities trending to zero should be disregarded. For example, the 
velocities at the lower end of the X=200.8 and 301 .6mm profiles decreased rapidly. One 
would expect the flow to be accelerated in this area by the high velocity plume shear 
layer. But the PIV did not capture these velocities due to low seeding or the delta-Ts 
being too large. 

Figure 62 plots entrained flowfield velocities at constant axial stations for the 
180° plane. The lines are not as smooth as the 0° plane. However, the characteristics just 
discussed about the flow in the 0° plane were all present in the 1 80° plane. Figure 63 
compares axial profiles from both planes to show the similarities. 
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Figure 56. Axial and Radial Stations of Subsequent Velocity Profiles for the Bell Nozzle. 
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0 Degree Plane, 8% Uncertainty 
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Figure 57. PIV Measured Velocities at Five Radial Stations in the 0° Plane for the Bell 

Nozzle. 
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180 Degree Plane, 8% Uncertainty 
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Figure 58. PIV Measured Velocities at Five Radial Stations in the 180° Plane for the Bell 

Nozzle. 
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0 and 180 Degree Planes, 8% Uncertainty 



Figure 59. Comparing PIV Measured Velocities for the 0 and 180° Planes at Three 

Radial Stations for the Bell Nozzle. 
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0 Degree Plane, 8% Uncertainty 
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Figure 60. PIV Measured Velocities at Five Axial Stations in the 0° Plane for the Bell 

Nozzle. 
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Figure 61 . PIV Measured Velocities at Five Additional Axial Stations in the 0° Plane for 

the Bell Nozzle. 


- 88 - 


FinalReport8.2 


11/15/2005 


180 Degree Plane, 8% Uncertainty 
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Figure 62. PIV Measured Velocities at Six Axial Stations in the 180° Plane for the Bell 

Nozzle. 
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Figure 63. Comparing PIV Measured Velocities for the 0 and 180° Planes at Axial 

Stations for the Bell Nozzle. 
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3.4.3. U nsteady Layer 

During testing, a region of unsteady flow was observed at the top of windows E to 
H in the Bell nozzle PIV data. In the average velocity plots, the unsteady layer appears 
as a low velocity region because the randomness of the velocity components drives the 
average toward zero. This low velocity region is seen at the top of windows E to H and 
the middle of window K in Figure 49 (and the top of Figure 33 for the Stratford nozzle). 
Evidence of it is also visible in the Bell nozzle velocity profiles in Figure 60 through 
Figure 63. An inflection exists in the velocity profiles at Y=350 and -350mm. This 
inflection marks the inner edge of the unsteady layer. 

This unsteadiness was evident in every delta-T dataset for these windows. 

Window K captured the entire width of this unsteady region. Six frames of instantaneous 
velocity vectors for window K are presented in Figure 64. The electronic media 
associated with this report contains an animation of velocity vectors for all 300 frames 
from one dataset for window K. 

During testing the unsteady layer was observed to be regular in its location and 
characteristics. The regions both above and below the layer had very steady uniform 
velocity vectors inward toward the plume; see Figure 64. Another unsteady layer 
appeared to be at the bottom of the FI 80 - HI 80 windows as well, although it was not as 
clearly delineated. Note that the layer was at the same radius at every axial station. It 
was not pulled toward the centerline by the radial inflow as one might expect. 

This unsteady layer is present in the Stratford PIV data as well. In the Stratford 
PIV data the unsteady layer was at a radius of approximately 260mm in the 0° plane and 
was only detected sporadically in the 1 80° plane. 

No explanation was found for this unsteady layer. Vortex shedding was 
considered. However, there was no hardware immediately upstream of the measurement 
windows that would induce vortex shedding. There were some flanges well upstream of 
the nozzle. But fog based flow visualization around the flanges did not indicate there was 
any vortex shedding off the flanges 

No indication of the unsteady layer was found in the hot film anemometer data for 
either the Stratford or Bell nozzles. 
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Figure 64. Velocity Vectors for Six of the 300 Instantaneous Frames of PIV Data from 

One Dataset of K Window. 
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3.4.4. PIV Measurement Repeatability 

The test schedule did not allow repeating PIV measurements in any of the 
windows. However, the overlapped regions of neighboring windows are effectively five 
regions where PIV data was obtained multiple times. Therefore, the repeatability of the 
PIV measured velocities was assessed by comparing the velocities in these overlapped 
regions. The difference in time between acquisition of the datasets was as little as a 
minute (within the same NTF run), or as much as a full day (in another NTF run). The 
repeatability of the PIV measurements were affected by at least the following: the 
variation of seeding quality, the accuracy of the camera positioning and the consistency 
of the test conditions (nozzle inlet and ambient). 

Repeatability was assessed by using the data from the regions where four 
windows overlapped. Five overlapped regions were analyzed for both nozzles. The 
overlapped regions were the coincident comers of the following sets of windows (see 
Figure 25 and Figure 46); A-B-E-F, B-C-F-G, C-D-G-H, B180-C180-F180-G180, and 
C180-D180-G180-H180. 

For the Stratford nozzle, the repeatability was assessed by comparing the 1 00 
microsecond delta-T datasets in the inner windows (A through D and B1 80 through 
D180) and 200 microsecond datasets in the outer windows (E through H and FI 80 
through HI 80). For the Bell nozzle, repeatability was assessed for the 0° plane by 
comparing the 200 microsecond delta-T datasets in the inner windows and the 500 
microsecond datasets for the outer windows. For the 180° plane 200 microsecond delta-T 
datasets were used for all windows. These delta-Ts were selected because they reflected 
the data selected by the merging technique discussed in Sections 3. 4. 1.2.2 and 3. 4. 1.2.2. 

Because the data points from the four windows of an overlapped region were not 
at the same locations in the global coordinate systems, the velocity data for each 
window’s dataset was linearly interpolated on to a common grid for that overlapped 
region. Therefore, each point in the common grid had four velocity values associated 
with it. A 95% confidence interval 5 was computed for this four sample set at each point. 
The confidence interval was then normalized by the average value at each point of the 
common grid. The minimum, maximum and average normalized confidence intervals for 
all the points in the common grids are presented in Table 7 through Table 10. The 
dimensional value of the average velocity components of the overlapped regions is also 
provided. 


Table 7. Stratford Nozzle U-Velocity Repeatability Summary 


Stratford Nozzle 

Normalized 95% Confidence Interval 

Average U Velocity 

Overlapping Windows 

MIN 

MAX 

Average 

(m/sec) 

A-B-E-F 

fFUiflSfj 

70 . 97 % 

28 . 65 % 

1.859 ± 0.53 

B-C-F-G 

HUSfl : 

26 . 57 % 

8 . 08 % 

2.566 ± 0.21 

C-D-G-H 

25 . 78 % 

60 . 46 % 

40 . 68 % 

2.807 ± 1.14 

B1 80-C1 80-F1 80-G1 80 

3 . 27 % 

19 . 89 % 

10 . 40 % 

2.084 ± 0.22 

C180-D1 80-G1 80-H1 80 

7 . 27 % 

38 . 44 % 

22 . 27 % 

2.554 ± 0.57 
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Table 8. Stratford Nozzle V-Velocity Repeatability Summary 


Stratford Nozzle 

Normalized 95% Confidence Interval 

Average V Velocity 

Overlapping Windows 

MIN 

MAX 

Average 

(m/sec) 

A-B-E-F 

43 . 23 % 

101 . 94 % 

62 . 17 % 

- 1.291 ± 0.80 

B-C-F-G 

23 . 77 % 

37 . 41 % 

31 . 16 % 

- 2.615 ± 0.82 

C-D-G-H 

4 . 14 % 

28 . 50 % 

19 . 77 % 

- 3.344 ± 0.66 

B180-C180-F180-G180 

15 . 95 % 

40 . 25 % 

23 . 64 % 

2.751 ± 0.65 

Cl 80-D1 80-G1 80-H1 80 

8 . 91 % 

39 . 80 % 

18 . 84 % 

3.462 ± 0.65 


Table 9. Bell Nozzle U-Velocity Repeatability Summary 


Bell Nozzle 

Normalized 95% Confidence Interval 

Average U Velocity 

Overlapping Windows 

MIN 

MAX 

Average 

(m/sec) 

A-B-E-F 

10 . 24 % 

12 . 58 % 

10 . 96 % 

2.121 ± 0.23 

B-C-F-G 

10 . 07 % 

16 . 00 % 

11 . 79 % 

2.058 ± 0.24 

C-D-G-H 

9 . 84 % 

12 . 01 % 

10 . 57 % 

2.279 ± 0.24 

B1 80-C1 80-F1 80-G1 80 

15 . 37 % 

66 . 92 % 

52 . 77 % 

1.433 ± 0.76 

C180-D1 80-G1 80-H1 80 

13 . 77 % 

36 . 57 % 

24 . 45 % 

1.912 ± 0.47 


Table 10. Bell Nozzle V-Velocity Repeatability Summary 


Bell Nozzle 

Normalized 95% Confidence Interval 

Average V Velocity 

Overlapping Windows 

MIN 

MAX 

Average 

(m/sec) 

A-B-E-F 



4 . 12 % 

- 2.730 ± 0.11 

B-C-F-G 



2 . 79 % 

- 3.634 ± 0.10 

C-D-G-H 

1 . 68 % 


2 . 41 % 

- 4.27 ± 0.10 

B1 80-C1 80-F1 80-G1 80 

5 . 81 % 

39 . 06 % 

21 . 00 % 

3.249 ± 0.67 

C180-D1 80-G1 80-H1 80 

1 . 60 % 

8 . 65 % 

4 . 39 % 

3.845 ± 0.17 


The Stratford nozzle PIV data was less repeatable than the Bell nozzle data. The 
Bell nozzle PIV data in the 0° plane was the most repeatable data. This is consistent with 
the general observations made elsewhere in this report. 
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3.4.5. Comparison of Stratford and Bell PIV Data 

In this section, observations made in previous sections from the PIV data plotted 
independently, are summarized first. Further observations are then made from concurrent 
plotting of the PIV data. 

3.4.5. 1 . Summary of Observations from PIV Data Plotted 
Independently 

The Stratford nozzle PIV data from the 1 80° plane was smoother than that from 
the 0° plane. The C window measurements were particularly poor due to seeding 
difficulties. The Bell nozzle PIV data from the 0° plane was smoother than that measured 
in the 180° plane. The Bell nozzle’s C window’s measurement quality was much 
improved relative to the Stratford’s C window, but is still of relatively lower quality than 
the neighboring windows within the Bell data. The highest PIV measured entrained flow 
velocity for the Stratford nozzle was in the 1 80° plane. For the Bell nozzle it was in the 
0° plane. 

Qualitatively the entrainment flowfields were similar. Both exhibit an 
acceleration of the entrained flow near the nozzle exit and a more gradual acceleration of 
the entire flowfield with axial station (as the flowfield approaches the plume capture 
pipe). 


In the regions near the nozzle exit, the maximum PIV measured entrained flow 
velocity for the Stratford nozzle was lower than the Bell nozzle. The Stratford PIV data 
indicates the U and V-velocity components of the entrained flow peaked about 25 and 
50mm downstream of the nozzle exit, respectively. The Bell PIV data indicates the U 
and V-velocity components peaked at about X=-10 and 10mm, respectively. The Bell 
nozzle flow separated just upstream of the nozzle exit inducing reverse flow into the 
nozzle. This reversed flow is the reason the Bell nozzle’s entrained flow velocities 
peaked closer to the nozzle exit than the Stratford. 

The effect of the acceleration near the nozzle exits on entrained flow is 
discemable in the PIV data in both the axial and radial directions. In the axial direction it 
is visible to about 125mm downstream of the nozzles. In the radial direction it is visible; 
for the Stratford, out to a radius of about 134mm, or 1 14mm beyond the nozzle outer 
diameter; for the Bell, out to a radius of about 227mm, or 145mm beyond the nozzle 
outer diameter. 


3. 4. 5. 2. Observations from PIV Data Plotted Concurrently 

The following paragraphs show that the two different nozzles induced different 
entrained flows near the nozzle exits, but the entrainment flowfield for the majority of 
area within the PIV measurement field was very similar. 

Figure 65 compares the Stratford 180° plane and the Bell 0° plane PIV data at four 
similar distances (‘dy’) from the respective nozzle outer diameters. The comparison was 
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done at these dy distances to compensate, approximately, for the different nozzle sizes. 
The Stratford 0° plane data is used in two places to compensate for bad or missing data 
from the 1 80° plane. 

The differences in the entrained flows near the nozzle exits are apparent between 
X=-100 and 125mm. The different axial locations of the peak measured velocities are 
clearly visible. Although not indicated by these curves, the Bell nozzle’s entrained flow 
velocities were, in fact, higher than the Stratford’s. These differences can be attributed to 
the following list of differences in expansion of nozzle flow by the Stratford and Bell 
nozzles. 

• The Stratford and Bell nozzles had inlet pressures of 10.35 and 19.65atm, 
respectively. 

• When expanded to ambient pressure (with isentropic equations), the plume Mach 
numbers were 2.2 and 2.6, respectively. 

• The Bell nozzle mass flow was approximately 2.3 times that of the Stratford 
nozzle. 

• The Bell nozzle’s plume cross-sectional and surface areas were 1.8 and 1.33 times 
that of the Stratford nozzle’s plume. 

• The Stratford nozzle was underexpanded, whereas, the Bell nozzle was 
overexpanded to the point of causing its flow to separate from the nozzle wall. 
This separation induced a sub ambient pressure region which induced a reverse 
flow into the Bell nozzle. 

Figure 65 also shows the similarities of the entrained flowfields for the majority 
of the area measured (i.e., outside of the region near the respective nozzle exits). At the 
second and third radii plotted (the red and green lines, respectively), the entrained flow 
velocities were very similar downstream of X= 125mm. By the third radii the entrained 
flow velocities were becoming more similar upstream of X= 125mm. These third radii 
are approaching the outer radii to which the Stratford and Bell nozzles affected the flow 
near the nozzle exit. At the fourth radii (the blue lines), well off the centerline and out of 
the nozzle exit effect, the entrained flow velocities were nearly identical. 

Figure 65, and this previous discussion of it, clearly shows that the Stratford and 
Bell nozzles did induce different entrained flows near the nozzle exit. But outside that 
region, the entrained flows were strikingly similar, especially, considering the long list of 
differences (given above) between the nozzles and their plumes. 

Now, in Figure 66, the entrained flow velocities for both nozzles are compared at 
constant axial stations. The Stratford PIV data is plotted on the right hand ordinate axis 
which was scaled to put both nozzles’ outer diameters at the same relative Y location. 

The peak velocities were different near the plume shear layer but the entrained flow 
velocities further from the centerline were quite similar (the differences at the outer edge 
of the profiles is discussed later in this subsection). The velocity profiles become more 
similar as they approach the plume capture pipe. 
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One exception was the X=300mm profile. Because the nozzle exits were not in 
the same axial location with respect to the plume capture pipe (the Stratford nozzle exit 
was 3 1 .75mm closer), the axial stations shown are not the same distance upstream from 
the plume capture pipe. To compensate for this for the 300mm profile, an additional 
profile was plotted from the Stratford nozzle data that was at the same location, relative 
to the plume capture pipe, as the Bell’s X=300mm profile. The agreement improved 
significantly. 

Figure 66 has again shown the entrained flow velocities, outside the plume shear 
layer, were strikingly similar. The velocity profiles became more similar as they 
approached the plume capture pipe. 

From these observations about the entrained flow’s velocities relative to the 
plume capture pipe, it appears the plume capture pipe sink was a fairly strong influence 
on the entrained flowfield. And that its effect on the entrained flow was similar for the 
two nozzles. 


Unsteady Layer Not From Vortex Shedding 

The unsteady layer’s effect on the PIV measured velocity magnitude is visible in 
the lower image of Figure 66. The inflection near the top of the velocity profiles marks 
the radial location of the unsteady layer’s lower edge. For the Stratford nozzle, in its 
global coordinate system, this was a radius of approximately 260mm; for the Bell nozzle, 
in its global coordinate system, approximately 340mm. Note, as indicated in Figure 66, 
this was at roughly the same radial distance relative to the nozzles’ outer diameters. This 
is further indication that the unsteadiness was not due to vortex shedding from some 
upstream hardware, but was due to some other, as of yet unexplained, flow instability. 
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Bell 0 and Stratford 180 Degree Planes, 8% Uncertainty 



— * — Bell, y=105.1mm, dy=23.8 

Strat, y=-44.5mm, dy=24.2 

— H — Bell, y=1 26.7mm, dy=45.4 
Strat, y=-65.9mm, dy=45.6 
— — Strat, y=66.0mm, dy=45.7 
— * — Bell, y=1 87.8mm, dy=106.5 
' Strat, y=-1 26.8mm, dy=106.5 

Strat, y=1 27.3mm, dy=107.0 

Bell, y=267.4mm, dy=186.1 
Strat, y=-206.6mm, dy=186.3 
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Figure 65. Comparison the Stratford and Bell Nozzle PIV Measured Velocities, for 
Constant Radial Stations, at Common Distances from the Nozzle Outer Diameter. 


- 98 - 


FinalReport8.2 


11/15/2005 
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Figure 66. Comparison the Stratford and Bell Nozzle PIV Measured Velocities, at 
Constant Axial Stations, at Common Distances from the Nozzle Outer Diameter. 
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3.4.6. PIV Data Files 

The electronic media associated with this report contains all PIV datasets and 
associated Excel files used to manipulate it. It also contains a few key TecPlot ‘layout’ 
files for plotting the data. The Bell and Stratford data structure and processing was very 
similar. The Bell nozzle files are described here. In these files the Stratford nozzle is at 
times referred to as the ‘Sonic’ nozzle. 

The PIV data for each nozzle is in its own folder. The majority of it was imported 
into Excel for processing, for either merging within a window or weaving an entire plane 
of data together. The top level Excel files ‘BellPIVDataSorted8.8pct.xls (and 
SonicPIVDataSorted8.8pct.xls) contain the merged and woven data for the 0°and 180° 
planes. It has worksheets containing the data and worksheets containing plots of the data. 
The data worksheets sort the PIV data at both constant axial and radial stations. This is 
the ‘woven’ data. The cells in those data worksheets reference another file, 
BellPIVData5.8pct.xls (SonicPIVData5.8pct.xls). It is in this file where the merging of 
the data in the inner rows of windows is preformed. The uncertainty criterion is set in 
cell U8 in the BellPIVData5.8pct.xls file. 

The worksheets containing plots in BellPIVDataSorted8.8pct.xls were used to create the 
line plots used within this report. They will be updated automatically if the uncertainty 
criterion is changed. 

The data worksheets of BellPIVDataSorted8.8pct were exported to create the TecPlot 
‘.dat’ data files. These .dat data files are included as well. 
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3.5. Hot Film Anemometer V elocity Data 

The hot film anemometer (HFA) data was pursued as a second set of data on the 
flow entrainment. The intent of the HFA was to: 

• To provide a dataset to fall back on if the PIV measurements were not successful. 

• To assist in interpreting the PIV measurement. 

• Obtain data over a larger area than was possible with the PIV technique. The 
measurement area in each plane was larger and data in the 90° plane was also 
obtained. 

The hot film anemometer wire was aligned to measure the radial component of 
velocity. No measurements were made with the HFA wire aligned to measure the 
horizontal, or U, component of velocity. The HFA technique only indicates a velocity 
magnitude. Direction of the velocity must be inferred from the probe orientation. The 
HFA velocity magnitude is presented as a positive value in all three planes. It is not 
assigned a direction (+ or -) based on the global coordinated systems Y axis. The HFA 
velocity data for the Stratford nozzle is discussed first. The HFA velocity data for the 
Bell nozzle is then discussed. The order and format in which the data are presented is 
very similar for both nozzles and is as follows. First, notes about the test data and the 
relative positions of the ‘windows’ are discussed. The HFA velocity field is then 
discussed. The HFA data for the two nozzle flowfields are compared. Then the HFA 
data is compared to the PIV V-velocity data. 

The HFA velocity data agreed well with the PIV data. The HFA data provided 
additional insight to the plume capture pipe sink’s effect on the entrained flow. 


3.5.1. Data Taken 

Position of Stratford Nozzle HFA Measurement Windows 

A schematic of the HFA measurement regions, or windows, for the Stratford 
nozzle is shown in Figure 67. Four windows were obtained in the 0°, 90° and 180° 
planes. 


Each window is 300 by 300mm. Measurements were taken every 20mm for a 
total of 256 measurement points in each window. Between runs the traverse was 
manually translated to the next window position, allowing for 1 00mm of overlap between 
the adjacent windows. 

The net HFA measurement area was much larger than the PIV measurement area. 
For comparison, Figure 69 indicates the PIV measurement area on a set of HFA data. 

The A window in each plane was positioned relative to the nozzle exit. The lower 
edge of the measurement region was set to be 25.4mm above the outer diameter of the 
nozzle at its exit. The left most edge of window A was limited by the Stratford nozzle 
shape. The left most edge of the A (and D) window was 40mm upstream of the nozzle 
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exit. The remaining windows in each plane were located by moving the traversing 
(and/or the rail system supporting it) a known distance relative to the A window. 

Coordinate transformation were performed to locate each windows in the 
Stratford nozzle’s global coordinate system. The transformations are noted in the Excel 
data files which contain the Stratford HFA test data. This Excel data files are described 
in Section 3.5.7. Note that the Stratford and Bell nozzle global coordinate systems are 
different. 


Position of Bell Nozzle HFA Measurement Windows 

A schematic of the HFA measurement regions, or windows, for the Bell nozzle is 
shown in Figure 68. Five windows were obtained in the 0° and 90° planes and four were 
obtained in the 180° plane. Each window is 300 by 300mm. Measurements were taken 
every 20mm for a total of 256 measurement points in each window. Between runs the 
traverse was manually translated to the next window position, allowing for 1 00mm of 
overlap between the adjacent windows. 

The net HFA measurement area was much larger than the PIV measurement area. 
For comparison, Figure 73 indicates the PIV measurement area on a set of HFA data. 

The A window in each plane was positioned relative to the nozzle exit. The lower 
edge of the measurement region was set to be 25.4mm above the outer diameter of the 
nozzle at its exit. The left most edge of the A (and D) window was 140mm upstream of 
the nozzle exit. The remaining windows in each plane were located by moving the 
traversing (and/or the rail system supporting it) a known distance relative to the A 
window. The right hand edge of the C window was 30mm from the plume capture pipe. 

Coordinate transformations were performed to locate each window in the Bell 
nozzle’s global coordinate system. The transformations are noted in the Excel data files 
which contain the Bell nozzle HFA test data. These Excel data files are described in 
Section 3.5.7. Note that the Stratford and Bell nozzle global coordinate systems are 
different. 

Notes from Run Logs 

Table 1 1 is a summary of the HFA test log. It was easier to change the nozzle test 
articles than it was to change the hot film traverse system and the rail system supporting 
it. Therefore, the test and instrumentation setup are listed in a chronological order instead 
of by nozzle. 
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Table 11. Summary of HFA Run Log 


Date 

Time 

Nozzle 

Window 

Plane 

Comment 

2/22/2005 

2:40 

Bell 

A 

0 

initial set-up with traverse on top of horizontal 
rail 

2/23/2005 

10:00 

Bell 

B 

0 

slide traverse, rail doesn't move 

2/23/2005 

1:50 

Bell 

C 

0 

slide traverse, rail doesn't move 

2/24/2005 

9:30 

Bell 

E 

0 

raise rail and slide traverse 

2/24/2005 

12:45 

Bell 

D 

0 

slide traverse 

2/25/2005 

10:20 

Bell 

A 

180 

lower rail and slide traverse 

2/25/2005 

12:30 

Bell 

B 

180 

slide traverse 

2/28/2005 

10:00 

Bell 

D 

180 

flip traverse to underside of rail to get planes 
close to floor 

2/28/2005 

1:30 

Bell 

E 

180 

slide traverse 

3/1/2005 

8:45 

Stratford 

A 

180 

install sonic nozzle - start with upside down 
traverse for all 180 

3/1/2005 

9:20 

Stratford 

B 

180 

OSC warning light at top 2 rows (Y=-45.7 and - 
65.7) in the 5 most downstream positions 
(X=380 to 460) 

3/1/2005 

1:00 

Stratford 

E 

180 

lower rail 

3/1/2005 

1:30 

Stratford 

D 

180 

slide traverse 

3/2/2005 

9:30 

Stratford 

A 

0 

probe wire broke while reconfiguring for 0 deg 
plane - calibrate new wire 

3/2/2005 

10:30 

Stratford 

B 

0 

slide traverse 

3/2/2005 

12:30 

Stratford 

E 

0 

raise rail 

3/2/2005 

1:00 

Stratford 

D 

0 

slide traverse 

3/3/2005 

8:30 

Stratford 

A 

90 

reconfigure for 90 deg plane - OSC warning 
light at (X=260, Y=45.7) position 

3/3/2005 

10:00 

Stratford 

B 

90 

OSC warning light at (X=440, Y=45.7) & 
(X=460, Y=45.7) positions 

3/3/2005 

11:00 

Stratford 

E 

90 

move rail 

3/3/2005 

1:00 

Stratford 

D 

90 

slide traverse 

3/4/2005 

10:30 

Bell 

A 

90 

install bell nozzle - facility data inadvertently 
not recorded for this run 

3/4/2005 

12:40 

Bell 

B 

90 

slide traverse 

3/4/2005 

2:15 

Bell 

C 

90 

slide traverse - OSC warning at (X=540, 
Y=106.6) 

3/7/2005 

9:00 

Bell 

E 

90 

move rail 

3/7/2005 

11:00 

Bell 

D 

90 

slide traverse 


The ‘OSC’ warning light mentioned in the comment column in Table 1 1 for 
several runs refers to a feature of the TSI Intelligent Flow Analyzer that indicated when 
the probe was in a highly oscillatory flow. It served as a warning to look at the turbulent 
intensity for those data points as well as real-time feedback to avoid damaging the probe. 
The warning light only went off when the probe was in or near the plume shear layer. 

The most significant item from the run logs was the probe change that occurred 
on 3/2/2005. The delicate wire of the hot film probe broke during a configuration change 
and had to be replaced. The replacement probe tip was calibrated in the same probe shaft 
and using the same calibration procedure. Testing proceeded without any noticeable 
problems. Unfortunately, differences between the probes later turned out to be 
significant. This only became apparent post test when comparing HFA and PIV data 
(Section 3.5.6). With nothing suspicious about the calibration of this second probe, its 
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disagreement with the rest of the data is attributed to the stack up of uncertainties 
associated with this HFA measurement technique. 



Figure 67. Schematic of HFA Datasets Acquired for the Stratford Nozzle. 
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Figure 68. Schematic of HFA Datasets Acquired for the Bell Nozzle. 
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3.5.2. Stratford Nozzle HFA Velocity Data 

Only one set of HFA measurements were made for each window. The HFA data 
from all windows was sorted in Excel to ‘weave’ the data together into one continuous 
set of data. The spacing of the HFA windows and data points within the windows was 
such that data points in the overlapped regions were at the same physical locations. For 
those coincident data points, an average value was calculated. 

Velocity Contour Plots 

Figure 69 presents the Stratford nozzle data for the 0 and 1 80° plane and Figure 
70 presents the same 0° plane of data but with the 90° plane of data. The HFA velocity 
magnitude contours appear quite similar in all three planes, particularly the 0 and 90° 
planes. The contours indicate the 1 80° plane had consistently lower velocities. As will 
be explained in Section 3.5.6, the velocities shown in the 180° plane are probably a more 
accurate representation of the entrained radial velocity component and the 0 and 90° 
plane HFA velocity magnitudes were probably higher than the true radial velocity. 

The indicated high velocity near the centerline is the plume shear layer. The 
indication of high velocity is correct, but the velocity magnitude indicated by the HFA 
probe is not accurate data. The probe was aligned to measure the radial component, but 
the axial component of the shear layer was two orders of magnitude larger than the radial 
component. See the vectors in Figure 34. With such large cross flows, the single wire 
HFA probe hardware used in this experiment could not provide accurate velocity data. 

The indicated high velocity on the right hand side of the measurements is the 
plume capture pipe sink. Figure 45 shows vectors of the sink flow near the pipe. The 
sink induced flow was predominately radial velocity flow in this location; therefore, these 
high indicated velocity magnitudes were not a result of axial cross flow. 

Velocity Line Plots 

Line plots of the HFA velocity magnitude at constant radial stations are now 
presented. The radial stations for these line plots are shown in Figure 71. The following 
line plots were created in an Excel file that is included in the electronic media associated 
with this report. The data plotted here can be extracted from that Excel file. This Excel 
file can be used to create additional line plots at alternate radial stations. Section 3.5.7 
discussed HFA Excel data files. 

Figure 72 plots the HFA velocity magnitude along several constant radial stations 
for the 0, 90 and 1 80° planes. The saw tooth patterns in the lines are the overlapping 
regions of adjacent windows. For the line plots, all values at overlapped points were 
plotted instead of an average value. As with the color contours above, the profiles in the 
0 and 90° planes were fairly similar. The profiles in the 1 80° plane had the same 
character as those in the 0 and 90 but the magnitude was consistently lower. Again, as 
explained in Section 3.5.6, the velocities shown in the 180° plane are probably a more 
accurate representation of the entrained radial velocity component and the 0 and 90° 
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plane HFA velocity magnitudes were probably higher than the true entrainment radial 
velocities. 

The Y=65.7 and 85.7mm profiles both show the acceleration near the nozzle exit. 
In the Y=65.7mm profiles, the higher indicated velocities to the right of X=250mm were 
a result of the shear layer cross flow and are suspect. The rest of the profiles show 
increasing velocity magnitude with axial station. For the outer profiles this acceleration 
was apparent over their entire length. This acceleration was mainly due to the sink 
created by the interaction of the plume and the plume capture pipe. The cause of the 
plume capture pipe sink is discussed in Section 3.4. 1.2. 3. 
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Figure 69. HFA Velocity Magnitude for Stratford Nozzle in the 0° and 180° Plane. 


- 107 - 



Y, Relative to Nozzle Centerline (mm) 


FinalReport8.2 


11/15/2005 


600 

500 

400 

300 

200 

100 

0 

100 

200 

300 

-400 

-500 

-600 


Hot Film Anemometer Data 
Velocity Magnitude 



I 1 1 I I I I I i ! i i I i 


Velocity Magnitude 

16 
14 
12 
10 
8 
6 
4 
2 
0 



1 I I I L 


0 200 400 600 800 1000 

X, Relative to Nozzle Exit (mm) 


1200 


Figure 70. HFA Velocity Magnitude for Stratford Nozzle in the 0° and 90° Plane. 
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Figure 71 . Radial Stations of HFA Velocity Profiles for the Stratford Nozzle. 
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Figure 72. HFA Velocity Magnitude at Constant Radial Stations for the Stratford 

Nozzle. 
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3.5.3. Bell Nozzle HFA Velocity Data 

Only one set of HFA measurements were made for each window. The HFA data from all 
windows was sorted in Excel to ‘weave’ the data together into one continuous set of data. 
The spacing of the HFA windows and data points within the windows was such that data 
points in the overlapped regions were at the same physical locations. For those 
coincident data points, an average value was calculated. 

Velocity Contour Plots 

Figure 73 presents the Bell nozzle data for the 0 and 180° plane and Figure 74 
presents the same 0° plane of data but with the 90° plane of data. The dark blue in the 
right hand comers in these figures are regions where no data was taken. The HFA 
velocity magnitude contours appear quite similar in all three planes. The contours 
indicate the 90° plane consistently had slightly higher velocity than the 0 and 1 80° planes. 
As will be explained in Section 3.5.6, the velocities shown in the 0 and 1 80° plane are 
probably a more accurate representation of the entrained radial velocity component and 
the 90° plane HFA velocity magnitudes were probably higher than the true radial 
velocity. 

In the Bell HFA data the plume shear layer was less of an influence than with the 
Stratford nozzle data. On the right hand side of data, near the centerline, a small portion 
of the plume shear layer was captured. The indication of high velocity in the shear layer 
is correct, but the velocity magnitude indicated by the HFA probe is not accurate data. 
The probe was aligned to measure the radial component, but the axial component of the 
shear layer in this region was significantly larger than the radial component. With large 
cross flows, the single wire HFA probe hardware used in this experiment could not 
provide accurate velocity data. 

The red bubble of high velocity on the right hand side of the measurements is the 
plume capture pipe sink. Figure 45 shows the vectors of the sink flow near the pipe. The 
sink induced flow was predominately radial velocity flow in this location, therefore, these 
high indicated velocity magnitudes were not a result of axial cross flow. 

Velocity Line Plots 

Line plots of the HFA velocity magnitude at constant radial stations are now 
presented. The radial stations for these line plots are shown in Figure 75. The following 
line plots were created in an Excel file that is included in the electronic media associated 
with this report. The data plotted here can be extracted from that Excel file. This Excel 
file can be used to create additional line plots at alternate radial stations. Section 3.5.7 
discussed HFA Excel data files. 

Figure 76 plots the HFA velocity magnitude along several constant radial stations 
for the 0, 90 and 1 80° planes. The saw tooth patterns in the lines are the overlapping 
regions of adjacent windows. For the line plots, all values at overlapped points were 
plotted instead of an average value. As with the color contours above, the profiles in the 
0 and 1 80° planes were fairly similar. The profiles in the 90° plane had the same 
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character as those in the 0 and 1 80 but the magnitude was consistently higher. Again, as 
explained in Section 3.5.6, the velocities shown in the 0 and 180° plane are probably a 
more accurate representation of the entrained radial velocity component and the 90° plane 
HFA velocity magnitudes are probably higher than the true entrainment radial velocity. 

The Y=106.7 and 126.7mm profiles both show the acceleration near the nozzle 
exit. In the Y=1 06.7mm profiles, the higher indicated velocities to the right of 
X=275mm were a result of the shear layer cross flow and are suspect. The right hand 
side of the Y=1 26.7mm profiles show rapidly increasing velocity. This is not a shear 
layer effect, but is the plume capture pipe induced radial flow. The rest of the profiles 
show increasing velocity magnitude with axial station. For the outer profiles this 
acceleration was apparent over their entire length. This acceleration was mainly due to 
the sink created by the interaction of the plume and the plume capture pipe. The cause of 
the plume capture pipe sink is discussed in Section 3. 4. 1.2. 3. 
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Figure 73. HFA Velocity Magnitude for Bell Nozzle in the 0° and 180° Plane. 
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Figure 74. HFA Velocity Magnitude for Bell Nozzle in the 0° and 90° Plane. 
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Figure 75. Radial Stations of HFA Velocity Profiles for the Bell Nozzle. 
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Figure 76. HFA Velocity Magnitude at Constant Radial Stations for the Bell Nozzle. 
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3.5.4. HFA Measurement Repeatability 

The test schedule did not allow repeating HFA measurements in any of the 
windows. However, the overlapped regions of neighboring windows are effectively five 
regions where HFA data was obtained multiple times. Therefore, the repeatability of the 
HFA measured velocities was assessed by comparing the velocities in these overlapped 
regions of the neighboring HFA datasets. Each window of data was acquired in a 
separate run of the NTF. The difference in time between acquisition of the datasets was 
as little as 30 minutes, or as much as a several days. The repeatability of the HFA 
measurements were affected by at least the consistency of the test conditions (nozzle inlet 
and ambient) and positional accuracy of the HFA hardware. 

Repeatability was assessed by using the data from the regions where three or four 
windows overlapped. Three overlapped regions were analyzed for Stratford nozzle and 
five for the Bell nozzle. The Stratford nozzle overlapped regions were the coincident 
comers of the following sets of windows (see Figure 67); A-B-D-E, A90-B90-D90-E90, 
and A180-B180-D180-E180. The Bell nozzle overlapped regions were the coincident 
comers of the following sets of windows (see Figure 68); A-B-D-E, B-C-E, A90-B90- 
D90-E90, B90-C90-E90, and A180-B180-D180-E180. 

Because the data points from the four windows of an overlapped region were at 
the same locations in the global coordinate systems, the velocity data for each window 
was compared with the positionally identical values of the other windows. Therefore, 
each data point location in the overlapping portion had three or four velocity values 
associated with it respectively. A 95% confidence interval 5 was computed for these three 
or four sample sets at each point, depending on the number of overlapping windows at 
each point. The confidence interval was then normalized by the average value at each 
point. The minimum, maximum and average normalized confidence intervals for all the 
points in the overlapping regions are presented in Table 12 and Table 13. The 
dimensional value of the average velocity magnitude of the overlapped regions is also 
provided. 


Table 12. Stratford Nozzle HFA Velocity Magnitude Repeatability Summary 


Stratford Nozzle 
Overlapping Windows 

Normalize 

MIN 

d 95% Confide 
MAX 

mce Interval 
Average 

Average Velocity 
Magnitude 
(m/sec) 

A-B-D-E 

A90-B90-D90-E90 
A180-B180-D1 80-El 80 


43 . 60 % 

106 . 35 % 

72 . 28 % 

19 . 92 % 

25 . 66 % 

36 . 10 % 

- 3.774 ± 0.75 
- 3.905 ± 1.00 
- 2.928 ± 1.06 


- 116 - 











FinalReport8.2 


11/15/2005 


Table 13. Bell Nozzle HFA Velocity Magnitude Repeatability Summary 






Average Velocity 

Bell Nozzle 

Normalized 95% Confidence Interval 

Magnitude 

Overlapping Windows 

MIN 

MAX 

Average 

(m/sec) 

A-B-D-E 

21 . 08 % 

66 . 40 % 

36 . 32 % 

- 2.566 ± 0.93 

B-C-E 

6 . 61 % 

60 . 61 % 

34 . 84 % 

- 3.821 ± 1.33 

A90-B90-D90-E90 

6 . 14 % 

42 . 07 % 

24 . 50 % 

- 3.179 ± 0.78 

B90-C90-E90 

7 . 87 % 

79 . 09 % 

26 . 69 % 

- 4.145 ± 1.11 

A1 80-B1 80-D1 80-El 80 

3 . 47 % 

91 . 36 % 

37 . 23 % 

- 2.489 ± 0.93 
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3.5.5. Comparison of Stratford and Bell HFA Data 

Although not shown here, the HFA data showed the acceleration of the entrained 
flow near the nozzle exits were different for the Stratford and Bell nozzles. These 
differences were similar to the differences shown by the PIV data (Section 3.4.5) for the 
near nozzle exit region of the Stratford and Bell nozzles. 

Figure 77 compares the outer radii profiles of the HFA velocities for the two 
nozzles. These outer profiles are quite similar in all three planes. The differences due to 
the acceleration of the entrained flow near the nozzle exits cannot be seen in these outer 
profiles. Both flowfields show the effect of the plume capture pipe sink; the acceleration 
of the entire field toward the plume capture pipe. 
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Figure 77. Comparison of Stratford and Bell Nozzle HFA Velocity Magnitude at the 

Outer Profiles. 
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3.5.6. Comparison of HFA Velocity Magnitude and PIV V-Velocity 

The HFA and PIV data are compared in the 0 and 180° planes for the Stratford 
and Bell nozzles. PIV data was not acquired on the 90° planes. Figure 78 plots the HFA 
velocity magnitude and PIV V-velocity for the Stratford nozzle 0° plane. The trends 
agree well but the HFA velocity is significantly higher than the PIV V-velocity. This is 
the only plane in which the two techniques do not agree. The other three planes agree 
exceptionally well (see Figure 79 through Figure 81). 

Two sets of HFA hardware (probe and wire combinations) were used to obtain the 
HFA data. The particular set used on the Stratford 0° plane was also used to acquire the 
Stratford nozzle’s 90° plane and the Bell nozzle’s 90° plane. The velocities recorded with 
this set of HFA hardware were consistently higher than those recorded (with the other set 
of HFA hardware) in the other three planes; the Stratford nozzle’s 180° and the Bell 
nozzle’s 0 and 180° planes. This is evident in the color contours; Figure 69 and Figure 
70, and Figure 73 and Figure 74; and the line plots, Figure 72 and Figure 76. 

It appears that HFA velocity magnitudes for the Stratford 0 and 90° plane and the 
Bell’s 90° plane are higher than the true entrainment radial velocities for these planes. 

This supposition is made because the HFA data in the Stratford nozzle’s 180° plane 
(Figure 79) and the Bell nozzle’s 0 and 180° planes (Figure 80 and Figure 81), all made 
with the other set of HFA hardware, were in excellent agreement with the PIV data. 

In these three planes, the only difference between results of the HFA and PIV 
techniques was at the outer most radial profiles of the Bell nozzle, where the PIV 
measurements captured the unsteady layer (see Section 3.4.3.) and the HFA did not. One 
hypothesis for this is that the frequency of the unsteadiness was relatively low. 

Sufficiently low such that the acquisition rate of the HFA technique, 1024 samples per 
second, was too high to capture the unsteadiness. The PIV technique acquired data at 
sufficiently slow rate, 1 5 samples per second, to capture some of the unsteadiness. 
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Stratford Nozzle 0 Degree Plane, PIV and HFA 
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Figure 78. Comparison of HFA Velocity Magnitude to PIV V-Velocity in the Stratford 

Nozzle 0° Plane. 
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Stratford Nozzle 180 Degree Plane, PIV and HFA 
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Figure 79. Comparison of HFA Velocity Magnitude to PIV V-Velocity in the Stratford 

Nozzle 180° Plane. 
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Bell Nozzle 0 Degree Plane, PIV and HFA 
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Figure 80. Comparison of HFA Velocity Magnitude to PIV V-Velocity in the Bell 

Nozzle 0° Plane. 


- 123 - 


HFA, Velocity Magnitude (m/s) HFA, Velocity Magnitude (m/s) 


Velocity (m/s) Velocity (m/s) 


FinalReport8.2 


11/15/2005 


Bell Nozzle 180 Degree Plane, PIV and HFA 
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Figure 81. Comparison of HFA Velocity Magnitude to PIV V-Velocity in the Bell 

Nozzle 180° Plane. 
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3.5.7. Entrained Flow Symmetry 

The PIV data for both nozzles (the Stratford, Figure 38; the Bell, Figure 59) 
showed that the entrained flowfield in the 0° planes was fairly similar to that in the 
respective 180° planes. One of the objectives of the HFA data was to determine if the 
entrained flowfields in the 90° planes was similar to that of 0 and 180° planes. That is, 
were the entrained flowfield symmetric about the nozzles? 

The HFA data showed that the character of the entrained flow in the 90° planes 
was similar to that in the 0 and 180° planes (Figure 72 and Figure 76). But as discussed 
in the section above, the HFA indicated velocity in the 90° planes was higher than that in 
the 0 and 180° planes. Both nozzles’ 90° planes were measured with the HFA hardware 
that appeared to record velocity magnitudes higher than the true velocity (see Section 
3.5.6). If in fact the 90° plane HFA data read high, and the real velocities were lower, the 
curves for the 90° planes (and the Stratford’s 0° plane) in Figure 72 and Figure 76 would 
move downward toward the curves for the other three planes. This would likely result in 
entrained flow that was more symmetric than that implied by Figure 72 and Figure 76. 

That said, the reader is reminded of the vortex at the inlet of the plume capture 
pipe in the 180° plane and the vortices attached to the west face of the test cabin (see 
Section 3.2). Both of which eliminate the possibility of truly axisymmetric entrained 
flowfields. 


3.5.8. HFA Data Files 

The electronic media associated with this report contains all HFA datasets and 
associated Excel files used to manipulate it. It also contains a few key TecPlot ‘layout’ 
files for plotting the data. The Bell and Stratford data structure and processing was very 
similar. The Bell nozzle files are described here. In these files the Stratford nozzle is at 
times referred to as the ‘Sonic’ nozzle. The HFA data for each nozzle is in its own 
folder. 


All of the raw data for each window was processed with two calibration curves, 
one for lower velocities and one for higher velocities. The velocities resulting from the 
‘low’ and ‘high’ calibrations are in neighboring columns in ‘Bellnearfield.xls’ 
(Stratnearfield.xls). Each widow’s data is contained in a separate worksheet. 

‘BellHFAdata6.xls’ (StratHFAdata2.xls) references Bellnearfield.xls and also has 
one worksheet for each window of data. In each worksheet a transition velocity (labeled 
‘filter’) is specified which sets the value below/above which the low/high calibrations 
were used. Most users of the data should have no reason to modify this transition value. 
‘BellHFAdata6.xls’ has worksheets (one per plane) that combine the data from all the 
windows in a plane. This worksheet is used to create the line plots presented in this 
report. 


‘BellHFAforTecPlot.xls’ (StratHFAforTecPlot.xls) references the window data in 
BellHFAdata6.xls. BellHFAforTecPlot.xls combines, or weaves, the data into one 
continuous data field for each plane. These woven data worksheets were exported to 
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create the TecPlot ‘.dat’ data files. The .dat data files are included in the electronic 
media. 

4. CONCLUSIONS 

The entrained flowfield velocities for the Bell and Stratford nozzles were captured 
with both particle image velocimetry (PIV) and hot film anemometry (HFA). The data 
acquired with the two different techniques exhibit excellent agreement in three of the four 
planes. For the one plane in which they did not agree as well, there is a reasonable 
explanation for the differences. 

In this test configuration there were two phenomena inducing the entrained flow. 
The first was the plume expansion at the nozzle exit. The second was the plume capture 
pipe sink. The data clearly shows that the acceleration of the entrained flows near the 
nozzle exits was different between the nozzles. The data also clearly shows that for the 
majority of the area in which the entrained flowfields were measured, the entrained flows 
were quite similar for both nozzles. This similarity was due to the sink that existed in the 
plume capture pipe. This sink resulted from the interaction of the plumes and the plume 
capture pipe inner wall. The sink induced flow from all directions toward the plume 
capture pipe inlet. 

The Stratford nozzle PIV data quality is acceptable. The PIV data acquisition 
techniques and test boundaries were being improved as the Stratford data was being 
acquired. Therefore, the Bell nozzle’s PIV data is the better quality data of the two PIV 
data sets. 

The PIV data indicated there was a region of unsteadiness in the entrained flow 
above the nozzle test articles. The source of the unsteady flow was not identified. 

The PIV and HFA data indicate the entrained flowfields in the 0, 90 and 180° 
planes were fairly similar. That is, the flow was generally symmetric about the nozzles. 
However, there were at least two large sources of asymmetry; the vortices and streaks of 
entrained flow from the base upstream of the nozzles and the vortex attached to the floor, 
going into the bottom of the plume capture pipe. 

Although the plume capture pipe sink had a significant effect on the entrained 
flowfield, the local effects of the nozzle flows could be discerned. Therefore, the data 
reported here should prove useful for benchmarking CFD codes for cold flow nozzle flow 
entrainment. 
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5. RECOMMENDATIONS AND LESSONS LEARNED 

In this section recommendations and lessons learned are presented in the general 
order that material was discussed in this report. 

Recommendation for CFD analysis of the test configuration: 

• The computational domain should be three dimensional. The entrained flow did 
contain three dimensional features. 

• However, the entrained flow was sufficiently symmetric that an axisymmetric 
computational domain is a plausible first approximation of the flow. 

• The plume capture pipe should be included in the computational domain. 
Attention should be paid to proper modeling of the plume capture pipe sink. 

Lessons Learned, general: 

• Parametric pretest CFD of the proposed test configurations was helpful in setting 
requirements for test configurations geometry and also for understanding the 
range of velocity magnitudes to be measured. 

• Data acquisition computers and equipment should be shielded from the high 
acoustic energy the plumes create. The ~ 125 to 133dB vibrated disk drives 
enough to make the computers inoperative. 

• The flow visualization with the fog, prior to obtaining any data, was very useful in 
understanding the large global flow features. 

• The plume is not fully captured by the pipe at low nozzle inlet pressures. For 
these experiments, the plume was fully captured when nozzle inlet pressures was 
greater than ~5atm. 

Lesson Learned related to the plume capture pipe sink effect on entrained flow: 

• Running two different nozzle configurations was very helpful in understanding 
where the entrainment effects due to the nozzle exit effect were dominate and 
where plume capture pipe sink effect was dominate. If only one nozzle had been 
run, then it would have been difficult to delineate the two effects. 

Recommendations related to the plume capture pipe sink effect on entrained flow: 

To reduce or eliminate the plume capture pipe sink effect: 

• Remove the plume capture pipe entirely, or 

• Remove more of it to increase distance between nozzle test article and pipe 
inlet, or 

• Perforate the pipe, with 25 or 50mm holes, downstream of its inlet to reduce 
the sink’s strength at the pipe inlet 

To more easily discern the effect of the plume capture pipe sink, obtain entrained 
flow velocities data: 

• With different distances between the nozzle exit and plume capture pipe inlet. 

• Without nozzles running, run the ejectors, currently in the pipe, at the same 
mass flow that the nozzles exhausted into the pipe. 
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Recommendations to better understand the effect of the Bell nozzle flow separation on 

entrained flow: 

• Obtain entrained flow velocity data for the current Bell nozzle with significantly 
lower nozzle inlet pressure, ~10.4atm. 

• Run a CFD analysis at the tested nozzle inlet pressure, then at a high enough 
nozzle inlet pressure to eliminate the nozzle flow separation. 

Recommendation to eliminate Bell nozzle flow separation: 

• Design a new test article with a lower area ratio. 

Recommendation to better understand the different effects on flow entrainment of the 

Stratford and Bell nozzles: 

• Run the Stratford at the same nozzle inlet pressure and/or mass flow as the Bell 
nozzle. 

Lessons Learned relative to entrainment flow data: 

• Acquiring both PIV and HFA data was useful. The PIV provided high density 
measurements. The HFA provided lower density measurement, but over a much 
larger area. This helped interpret the ‘big picture’ of the flow physics and helped 
in understanding the PIV data better. 

Lessons Learned relative to the PIV data acquisition: 

• Have the openings in the building far from the test article and measurement 
region. Openings near the measurement fields induced additional fluctuations in 
the entrained flow. 

• Place the seeder far from the PIV measurement field so as to seed the entrained 
air as evenly as possible. For this test, the seeder by the garage door worked well. 

• Perform all seeding parametrics and diagnostics at reduced nozzle inlet pressure 
to conserve facility air supply. 

• Two scheduled data acquisition periods, with several days between them, enabled 
problems discovered in the first period to be fixed without wasting valuable PIV 
test time. 

• Do not try to acquire the highest priority data in the first test runs. Use the first 
runs to understand the PIV system requirements for obtaining good data and to 
understand the structure of the nozzle flow. 

Recommendations for future PIV data acquisition: 

• The CCD camera used for PIV testing was a 2 Mega Pixel camera. Use of new 
generation, more sensitive with a higher number of active pixels, CCD cameras 
such as 4 or 1 1 Mega Pixel cameras would allow for the capture of a much larger 
field of view and negate the need for merging smaller FOV PIV data sets. 

• The need to illuminate a large field of view would also require more powerful 
dual cavity YAG lasers than the 120mJ per pulse one used here. 

• Use multiple PIV cameras side by side to double the effective field of view, and, 
therefore, data acquired. 

• Shield the cameras and computers from the acoustic environment generated by 
the plumes. 
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• Use multiple seeders to ensure optimal seed particle uniformity and density. 

Recommendations to improve the HFA data acquisition: 

• Obtain and use a two component instrument so both U and V- velocity 
components could be captured at the same time. 

• Have multiple probes on hand prior to starting the test. They break easily. 

• Have probes calibrated prior to testing. 

• Have real time feed back during calibration of HFA probes. 

Lesson Learned relative to capturing flow visualization with the digital video camera and 
low power laser: 

• Insufficient light energy was available if the camera was position to capture ‘back 
scattered’ light. That is, if the laser and camera were on the same side of the 
feature being illuminated with light scattered backward toward the camera. 
Sufficient light energy was available when the camera was positioned to capture 
‘forward scattered’ light. That is, if the camera was on the other side of the 
feature being illuminated with light scattered forward toward the camera. More 
light is scattered forward than back. 
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7. APPENDICIES 

Appendix 1. Nozzle Test Article Geometries and Bell Nozzle Wall Pressure Data. 

The nozzle internal geometries are presented in Tables A1 and A2. The Stratford 
nozzle internal geometry is the as-designed contour for a throat of 35.56mm. The Bell 
nozzle internal geometry was extracted from an Initial Graphics Exchange Specification 
file, refined in a grid generator and printed out. Bell nozzle wall static pressure 
measurement locations are documented in Table A3. 
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3.5.7. Entrained Flow Symmetry 

The PIV data for both nozzles (the Stratford, Figure 38; the Bell, Figure 59) 
showed that the entrained flowfield in the 0° planes was fairly similar to that in the 
respective 180° planes. One of the objectives of the HFA data was to determine if the 
entrained flowfields in the 90° planes was similar to that of 0 and 180° planes. That is, 
were the entrained flowfield symmetric about the nozzles? 

The HFA data showed that the character of the entrained flow in the 90° planes 
was similar to that in the 0 and 180° planes (Figure 72 and Figure 76). But as discussed 
in the section above, the HFA indicated velocity in the 90° planes was higher than that in 
the 0 and 180° planes. Both nozzles’ 90° planes were measured with the HFA hardware 
that appeared to record velocity magnitudes higher than the true velocity (see Section 
3.5.6). If in fact the 90° plane HFA data read high, and the real velocities were lower, the 
curves for the 90° planes (and the Stratford’s 0° plane) in Figure 72 and Figure 76 would 
move downward toward the curves for the other three planes. This would likely result in 
entrained flow that was more symmetric than that implied by Figure 72 and Figure 76. 

That said, the reader is reminded of the vortex at the inlet of the plume capture 
pipe in the 180° plane and the vortices attached to the west face of the test cabin (see 
Section 3.2). Both of which eliminate the possibility of truly axisymmetric entrained 
flowfields. 


3.5.8. HFA Data Files 

The electronic media associated with this report contains all HFA datasets and 
associated Excel files used to manipulate it. It also contains a few key TecPlot ‘layout’ 
files for plotting the data. The Bell and Stratford data structure and processing was very 
similar. The Bell nozzle files are described here. In these files the Stratford nozzle is at 
times referred to as the ‘Sonic’ nozzle. The HFA data for each nozzle is in its own 
folder. 


All of the raw data for each window was processed with two calibration curves, 
one for lower velocities and one for higher velocities. The velocities resulting from the 
‘low’ and ‘high’ calibrations are in neighboring columns in ‘Bellnearfield.xls’ 
(Stratnearfield.xls). Each widow’s data is contained in a separate worksheet. 

‘BellHFAdata6.xls’ (StratHFAdata2.xls) references Bellnearfield.xls and also has 
one worksheet for each window of data. In each worksheet a transition velocity (labeled 
‘filter’) is specified which sets the value below/above which the low/high calibrations 
were used. Most users of the data should have no reason to modify this transition value. 
‘BellHFAdata6.xls’ has worksheets (one per plane) that combine the data from all the 
windows in a plane. This worksheet is used to create the line plots presented in this 
report. 


‘BellHFAforTecPlot.xls’ (StratHFAforTecPlot.xls) references the window data in 
BellHFAdata6.xls. BellHFAforTecPlot.xls combines, or weaves, the data into one 
continuous data field for each plane. These woven data worksheets were exported to 
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create the TecPlot ‘.dat’ data files. The .dat data files are included in the electronic 
media. 

4. CONCLUSIONS 

The entrained flowfield velocities for the Bell and Stratford nozzles were captured 
with both particle image velocimetry (PIV) and hot film anemometry (HFA). The data 
acquired with the two different techniques exhibit excellent agreement in three of the four 
planes. For the one plane in which they did not agree as well, there is a reasonable 
explanation for the differences. 

In this test configuration there were two phenomena inducing the entrained flow. 
The first was the plume expansion at the nozzle exit. The second was the plume capture 
pipe sink. The data clearly shows that the acceleration of the entrained flows near the 
nozzle exits was different between the nozzles. The data also clearly shows that for the 
majority of the area in which the entrained flowfields were measured, the entrained flows 
were quite similar for both nozzles. This similarity was due to the sink that existed in the 
plume capture pipe. This sink resulted from the interaction of the plumes and the plume 
capture pipe inner wall. The sink induced flow from all directions toward the plume 
capture pipe inlet. 

The Stratford nozzle PIV data quality is acceptable. The PIV data acquisition 
techniques and test boundaries were being improved as the Stratford data was being 
acquired. Therefore, the Bell nozzle’s PIV data is the better quality data of the two PIV 
data sets. 

The PIV data indicated there was a region of unsteadiness in the entrained flow 
above the nozzle test articles. The source of the unsteady flow was not identified. 

The PIV and HFA data indicate the entrained flowfields in the 0, 90 and 180° 
planes were fairly similar. That is, the flow was generally symmetric about the nozzles. 
However, there were at least two large sources of asymmetry; the vortices and streaks of 
entrained flow from the base upstream of the nozzles and the vortex attached to the floor, 
going into the bottom of the plume capture pipe. 

Although the plume capture pipe sink had a significant effect on the entrained 
flowfield, the local effects of the nozzle flows could be discerned. Therefore, the data 
reported here should prove useful for benchmarking CFD codes for cold flow nozzle flow 
entrainment. 
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5. RECOMMENDATIONS AND LESSONS LEARNED 

In this section recommendations and lessons learned are presented in the general 
order that material was discussed in this report. 

Recommendation for CFD analysis of the test configuration: 

• The computational domain should be three dimensional. The entrained flow did 
contain three dimensional features. 

• However, the entrained flow was sufficiently symmetric that an axisymmetric 
computational domain is a plausible first approximation of the flow. 

• The plume capture pipe should be included in the computational domain. 
Attention should be paid to proper modeling of the plume capture pipe sink. 

Lessons Learned, general: 

• Parametric pretest CFD of the proposed test configurations was helpful in setting 
requirements for test configurations geometry and also for understanding the 
range of velocity magnitudes to be measured. 

• Data acquisition computers and equipment should be shielded from the high 
acoustic energy the plumes create. The -125 to 133dB vibrated disk drives 
enough to make the computers inoperative. 

• The flow visualization with the fog, prior to obtaining any data, was very useful in 
understanding the large global flow features. 

• The plume is not fully captured by the pipe at low nozzle inlet pressures. For 
these experiments, the plume was fully captured when nozzle inlet pressures was 
greater than ~5atm. 

Lesson Learned related to the plume capture pipe sink effect on entrained flow: 

• Running two different nozzle configurations was very helpful in understanding 
where the entrainment effects due to the nozzle exit effect were dominate and 
where plume capture pipe sink effect was dominate. If only one nozzle had been 
run, then it would have been difficult to delineate the two effects. 

Recommendations related to the plume capture pipe sink effect on entrained flow: 

To reduce or eliminate the plume capture pipe sink effect: 

• Remove the plume capture pipe entirely, or 

• Remove more of it to increase distance between nozzle test article and pipe 
inlet, or 

• Perforate the pipe, with 25 or 50mm holes, downstream of its inlet to reduce 
the sink’s strength at the pipe inlet 

To more easily discern the effect of the plume capture pipe sink, obtain entrained 
flow velocities data: 

• With different distances between the nozzle exit and plume capture pipe inlet. 

• Without nozzles running, run the ejectors, currently in the pipe, at the same 
mass flow that the nozzles exhausted into the pipe. 
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Recommendations to better understand the effect of the Bell nozzle flow separation on 

entrained flow: 

• Obtain entrained flow velocity data for the current Bell nozzle with significantly 
lower nozzle inlet pressure, ~10.4atm. 

• Run a CFD analysis at the tested nozzle inlet pressure, then at a high enough 
nozzle inlet pressure to eliminate the nozzle flow separation. 

Recommendation to eliminate Bell nozzle flow separation: 

• Design a new test article with a lower area ratio. 

Recommendation to better understand the different effects on flow entrainment of the 

Stratford and Bell nozzles: 

• Run the Stratford at the same nozzle inlet pressure and/or mass flow as the Bell 
nozzle. 

Lessons Learned relative to entrainment flow data: 

• Acquiring both PIV and HFA data was useful. The PIV provided high density 
measurements. The HFA provided lower density measurement, but over a much 
larger area. This helped interpret the ‘big picture’ of the flow physics and helped 
in understanding the PIV data better. 

Lessons Learned relative to the PIV data acquisition: 

• Have the openings in the building far from the test article and measurement 
region. Openings near the measurement fields induced additional fluctuations in 
the entrained flow. 

• Place the seeder far from the PIV measurement field so as to seed the entrained 
air as evenly as possible. For this test, the seeder by the garage door worked well. 

• Perform all seeding parametrics and diagnostics at reduced nozzle inlet pressure 
to conserve facility air supply. 

• Two scheduled data acquisition periods, with several days between them, enabled 
problems discovered in the first period to be fixed without wasting valuable PIV 
test time. 

• Do not try to acquire the highest priority data in the first test runs. Use the first 
runs to understand the PIV system requirements for obtaining good data and to 
understand the structure of the nozzle flow. 

Recommendations for future PIV data acquisition: 

• The CCD camera used for PIV testing was a 2 Mega Pixel camera. Use of new 
generation, more sensitive with a higher number of active pixels, CCD cameras 
such as 4 or 1 1 Mega Pixel cameras would allow for the capture of a much larger 
field of view and negate the need for merging smaller FOV PIV data sets. 

• The need to illuminate a large field of view would also require more powerful 
dual cavity YAG lasers than the 120mJ per pulse one used here. 

• Use multiple PIV cameras side by side to double the effective field of view, and, 
therefore, data acquired. 

• Shield the cameras and computers from the acoustic environment generated by 
the plumes. 
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• Use multiple seeders to ensure optimal seed particle uniformity and density. 

Recommendations to improve the HFA data acquisition: 

• Obtain and use a two component instrument so both U and V-velocity 
components could be captured at the same time. 

• Have multiple probes on hand prior to starting the test. They break easily. 

• Have probes calibrated prior to testing. 

• Have real time feed back during calibration of HFA probes. 

Lesson Learned relative to capturing flow visualization with the digital video camera and 
low power laser: 

• Insufficient light energy was available if the camera was position to capture ‘back 
scattered’ light. That is, if the laser and camera were on the same side of the 
feature being illuminated with light scattered backward toward the camera. 
Sufficient light energy was available when the camera was positioned to capture 
‘forward scattered’ light. That is, if the camera was on the other side of the 
feature being illuminated with light scattered forward toward the camera. More 
light is scattered forward than back. 
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7. APPENDICIES 

Appendix 1. Nozzle Test Article Geometries and Bell Nozzle Wall Pressure Data. 

The nozzle internal geometries are presented in Tables A1 and A2. The Stratford 
nozzle internal geometry is the as-designed contour for a throat of 35.56mm. The Bell 
nozzle internal geometry was extracted from an Initial Graphics Exchange Specification 
file, refined in a grid generator and printed out. Bell nozzle wall static pressure 
measurement locations are documented in Table A3. 
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Table A2. Bell No zzle Internal, As-Design ed, Wall Geometry. 
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Table A3. Bell Nozzle Wall Static Pressure Measurement Locations 


Wall Pstatic 
Measurement 
Label 

X(in) Relative 
to Nozzle 
Throat 

X(mm) 
Relative to 
Nozzle Throat 

X(mm) Relative 
to Global 
Coordinate 
System 

PS2-SR-1 

0.2931 

7.4447 

-138.5310 

PS2-SR-2 

0.5592 

14.2037 

-131.7721 

PS2-SR-3 

0.8868 

22.5247 

-123.4510 

PS2-SR-4 

1.176 

29.8704 

-116.1054 

PS2-SR-5 

1.5042 

38.2067 

-107.7691 

PS1-SR-6 

1.816 

46.1264 

-99.8494 

PS1-SR-7 

2.1121 

53.6473 

-92.3284 

PS1-SR-8 

2.4223 

61.5264 

-84.4493 

PS1-SR-9 

2.6853 

68.2066 

-77.7691 

PS1-SR-10 

3.005 

76.3270 

-69.6488 

PS1-SR-1 1 

3.3094 

84.0588 

-61.9170 

PS1-SR-12 

3.6134 

91.7804 

-54.1954 

PS1-SR-13 

3.9094 

99.2988 

-46.6770 

PS1-SR-14 

4.2259 

107.3379 

-38.6379 

PS1-SR-15 

4.5235 

114.8969 

-31.0789 

PS1-SR-16 

4.8306 

122.6972 

-23.2785 

PS1-SR-17 

5.1124 

129.8550 

-16.1208 

PS1-SR-18 

5.4132 

137.4953 

-8.4805 

PS2-SR-19 

0.2777 

7.0536 

-138.9222 

PS2-SR-20 

0.5532 

14.0513 

-131.9245 

PS2-SR-21 

0.8531 

21.6687 

-124.3070 

PS2-SR-22 

1.1539 

29.3091 

-116.6667 

PS2-SR-23 

1.4751 

37.4675 

-108.5082 

PS1-SR-24 

1.7759 

45.1079 

-100.8679 

PS1-SR-25 

2.0968 

53.2587 

-92.7170 

PS1-SR-26 

2.4101 

61.2165 

-84.7592 

PS1-SR-27 

2.714 

68.9356 

-77.0402 

PS1-SR-28 

2.9896 

75.9358 

-70.0399 

PS1-SR-29 

3.3171 

84.2543 

-61.7214 

PS1-SR-30 

3.6258 

92.0953 

-53.8804 

PS1-SR-31 

3.9045 

99.1743 

-46.8015 

PS1-SR-32 

4.2053 

106.8146 

-39.1611 

PS1-SR-33 

4.4966 

114.2136 

-31.7621 

PS1-SR-34 

4.8068 

122.0927 

-23.8830 

PS1-SR-35 

5.1155 

129.9337 

-16.0421 

PS1-SR-36 

5.4179 

137.6147 

-8.3611 

PS3-SR-37 

-1 .80 

-45.7200 

-191.6958 

PS3-SR-38 

-1.80 

-45.7200 

-191.6958 

PS3-SR-39 

-1.80 

-45.7200 

-191.6958 
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Table A4. Bell Nozzle Wall Pressures 





Nozzle Pressure Ratio 

Pressure 

Tap 

Axial Station 
Relative to 
Nozzle Exit 
(mm) 

Axial Station 
Relative to 
Nozzle 
Throat (mm) 

10.41 

19.83 

19.87 

19.89 

19.15 




Pwall/Pnozzlejnlet 

PS2-SR-1 

-138.5310 

7.4447 


0.0984 

0.0983 

0.0986 

0.0986 

PS2-SR-2 


14.2037 



■■Ifei 

0.0861 

0.0861 

PS2-SR-3 


22.5247 





0.0695 

PS2-SR-4 


29.8704 

0.0541 

0.0540 

0.0540 

■masi 

■iBlgHi 

PS2-SR-5 


38.2067 


■neitM 



l—w 

PS1-SR-6 

-99.8494 

46.1264 

■iliEliM 

■iiiwai 

■lillEM 



PS1-SR-7 

-92.3284 

53.6473 

IltBl 


■iliHM 


WSSEM 

PS1-SR-8 

-84.4493 

61.5264 

0.0898 

0.0285 

0.0285 

■lllfesa 

EES9 

PS1-SR-9 

-77.7691 

68.2066 


■UIWM 

0.0253 

■lliwsl 

liTi^cl 

PS1-SR-10 

-69.6488 

76.3270 

■item 





PS1-SR-11 

-61.9170 

84.0588 

■IIIKWIM 

■iliHiEB 

0.0203 

0.0203 


PS1-SR-12 

-54.1954 

91.7804 

■lliliVJeB 

■IIIIHSM 

0.0183 

0.0184 

■HllliM 

PS1-SR-13 

-46.6770 

99.2988 


■iiiiriM 

0.0171 


■iliiMB 

PS1-SR-14 

-38.6379 

107.3379 

0.0924 

0.0157 

0.0157 

0.0158 

0.0158 


-31.0789 

114.8969 

0.0926 

0.0169 

0.0164 

0.0160 

0.0340 

PS1-SR-16 

-23.2785 


0 0928 

0.0387 

0.0384 

0.0381 

0.0437 

PS1-SR-17 

-16.1208 

129.8550 

0.0928 

0.0469 


■TOM 

0.0490 

PS1-SR-18 

-8.4805 


0.0930 

ESI 


■iliW 

0.0472 

PS2-SR-19 

-138.9222 

7.0536 




■iiieniM 

0.0987 

PS2-SR-20 

-131.9245 


0.0855 


0.0844 

0.0845 

0.0843 

PS2-SR-21 

-124.3070 

21.6687 

0.0690 

0.0686 


0.0688 

0.0688 

PS2-SR-22 

-116.6667 

29.3091 

0.0550 

0.0550 



0.0551 

PS2-SR-23 

-108.5082 

37.4675 


0.0455 

0.0455 

0.0457 

0.0457 

PS1-SR-24 

-100.8679 

45.1079 

■ill HbM 

ESI 

0.0384 

0.0384 

eemi 

PS1-SR-25 

-92.7170 

53.2587 

■IIIEffil 

mmstm 



mBBsm 

PS1-SR-26 

-84.7592 

61.2165 


0.0279 

0.0279 

0.0280 

0.0280 

PS1-SR-27 

-77.0402 

68.9356 


0.0244 



0.0245 

PS1-SR-28 

-70.0399 

75.9358 

0.0916 

0.0224 

0.0224 

0.0225 

0.0225 

PS1-SR-29 

-61.7214 

84.2543 

0.0922 

0.0195 

0.0195 


0.0196 

PS1-SR-30 

-53.8804 

92.0953 

0.0924 

0.0185 

0.0185 

0.0185 

0.0185 

PS1-SR-31 

-46.8015 

99.1743 


0.0172 

■ihIm 

0.0172 

0.0172 

PS1-SR-32 

-39.1611 

106.8146 

MUM 

0.0158 

■iliiW 

0.0159 


PS1-SR-33 

-31.7621 

114.2136 

0.0927 

0.0175 

0 0164 

0.0159 


PS1-SR-34 

-23.8830 

122.0927 

0.0930 

0.0393 


0.0389 

0.0441 | 

PS1-SR-35 

-16.0421 

129.9337 

0.0932 

0.0449 

0.0449 

0.0448 

wesm 

PS1-SR-36 

-8.3611 

137.6147 



0.0470 

0.0469 


PS3-SR-37 

-191.6958 

-45.7200 

0.9773 

0.9761 

0.9753 

0.9753 

WBEM 

PS3-SR-38 

-191.6958 

-45.7200 



0.9769 

0.9767 

0.9766 | 

PS3-SR-39 

-191.6958 

-45.7200 






Nozzle Base 1A 

0.0000 

145.9758 

0.09453 



rrrwm 

0.049552 

0.051429 

Nozzle Base 2A 

^MTiTiTiTiW^ 

145.9758 

0.09486 

0.049802 

0.049717 

0.04966 

0.051536 

Nozzle Base IB 


145.9758 

0.09466 

0.049686 

0.049567 

0.049541 

0.051412 

Nozzle Base 2B 


145.9758 

0.09570 

0.050233 

0.050143 

0.050094 

0.052026 

Nozzle Base 1C 

0.0000 

145.9758 

0.09503 

0.049921 

0.049809 

0.049771 

0.051692 

Nozzle Base 2C 

0.0000 

145.9758 

0.09510 

0.049944 

0.04987 

0.049791 

0.051705 
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